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PART I 


REVIEW AND EVALUATION OF RECENT DEVELOPMENTS 
IN MERLICK INLET DYNAMIC FLOW DISTORTION 
PREDICTION 



SUMMARY 


A brief review of developments in the Melick method 
of inlet flow dynamic distortion prediction by statistical 
means is provided. These developments include the general 
Melick approach with full dynamic measurements, a limited 
dynamic measurement approach, and a turbulence modelling 
approach which requires no dynamic rms pressure fluctuation 
measurements. These modifications are briefly evaluated 
by comparing predicted and measured peak instantaneous 
distortion levels from provisional inlet data sets. 

A nonlinear mean-line fallowing vortex model is pro- 
posed and evaluated as a potential criterion for improving 
the peak instantaneous distortion map generated from 
the conventional linear vortex of the Melick method. The 
model is simplified to a series of linear vortex segments 
which lay along the mean line. Maps generated with this new 
approach are compared with conventionally generated maps , 
as well as measured peak instantaneous maps. 

Results of the developments and modifications discussed 
compare well with experimental measurements, both in the 
prediction of peak instantaneous distortion levels, and the 
peak instantaneous maps. Inlet data sets include subsonic , 
transonic, and supersonic inlets under various flight 
conditions. The methods discussed can be used in preliminary 
inlet design phases in the i nterest of reducing development 


costs . 
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INTRODUCTION 


I n let: turbulence and other flow nonuniformities have 
long been known to significantly affect the operational 
stability of gas turbine engines, especially in high per- 
formance military aircraft. This inlet flow distortion is 
traditionally measured at the compressor face of the engine 
with an array of total pressure probes mounted an rakes. 

The time-averaged steady-state pressures at each of the 
probe locations are processed and combined in such a way as 
to generate various steady-state distortion factors and an 
engine face pressure contour map (see Table 1 and Figure 1, 
respectively). These then correlate to engine surge margins . 

The distortion problem is intensified by the time var- 
iant component of the total distortion level. Random fluc- 
tuations in the total pressure measurements can generate 
instantaneous distortion levels which can induce engine 
surges even when the steady-state component is well below 
compressor stall margins. It becomes important, therefore, 
to be able to predict the most probable peak instantaneous 
(dynamic) distortion level early in the inlet design effort. 

One method of determining the dynamic distortion level 
of an inlet is to use an array of high response total pres- 
sure probes, with an extensive inventory of support instru- 
mentation and computationa 1 equipment to record time histo- 
ries of the pressure fluctuations for each of the probes. 
These data are then screened, using the Dynamic Data Editing 
and Computing (DYNADEC) system, to determine an experimental 
peak distortion level using the same definitions as the 
steady-state case. This method is generally quite accurate, 
compared to statistical methods described later, but it is 
also extremely expensive in terms of instrumentation and 
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computational requirements. In order to reduce the cost of 
inlet distortion tests, several statistical methods have 
been developed to predict the dynamic distortion component, 
given the steady-state distortion and limited dynamic data 
(ref. 1, 2, 3) . 

OF the many statistical methods of predicting dynamic 
distortion levels, the most efficient is Melick modelling 
approach. In the Melick method, it is postulated that the 
dynamic disturbances in the inlet flow can be modelled by 
the pressure disturbances resulting from a series of random- 
ly distributed vortices convecting through the inlet duct. 
Filtered and unfiltered root mean square (rms) total pres- 
sure fluctuation levels are used to identify the main varia- 
bles in this vortex flow model (ref. 2, 3, 4). 

The main advantages of the Melick method include low 
cost relative to other techniques, as well as the fact that 
it can be used online, while the test is in progress. It 
has been shown that further cost reduction can be attained 
by reducing the quantity of dynamic data (ref. 2, 3}„ In 

fact, Chen (ref. 3) has derived and demonstrated a new tech- 
nique for predicting the peak distortion levels with only 
the steady-state distortion data, that is, with no dynamic 
data . 

One of the main disadvantages of the Melick modelling 
approach is it is not as accurate as some methods in the 
generation of the peak dynamic distortion patterns in the 
engine face contour map. This is due primarily to a limita- 
tion in the vortex flow model, namely, the use of a single 
linear vortex in the generation of the peak instantaneous 
pressure array. More specifically, the peak instantaneous 
pressure array is computed by placing a linear vortex along 
a portion of the mean pressure line in such a way to amplify 
the distortion pattern pressures (fig. 3). The size and 
strength of this vortex Is determined as a function of the 
most probable peak instantaneous distortion level. It has 
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been suggested that; a new concept: in vortex modelling could 
improve the accuracy of the predicted peak dynamic distor- 
tion pattern (ref, 2). 

It is apparent that the mean pressure line in distor- 
tion patterns is not generally straight. In most cases, the 
mean line can be seen to arc across the engine face, fre- 
quently forming a distorted ring. One passible solution to 
the vortex modelling problem is ta replace the single 
straight vortex oriented along a portion of the mean line 
with a curved mean- 1 ine-fo 1 lowing vortex. This nonlinear 
vortex (or vortex ring, where applicable] could provide a 
more accurate amplification of the pressure levels in the 
vicinity of the vortex. 

In the present work, the concept of replacing the 
linear vortex model of the traditional Melick model with a 
nonlinear mean-line-following vortex is proposed and 
evaluated. For the purposes of demonstrating the concept, 
this new model is simplified by breaking the vortex into 
a series of vortex segments, one segment for each of the 
probe rakes (fig. 4] . The radius and strength of these 
vortex segments is retained from the original Melick dynamic 
data matching process. 

The results of the present method are compared to the 
original single linear vortex model, as well as the DYNADEC 
results, for a variety of data sets. These data include 
example subsonic, transonic, and supersonic inlet configu- 
rations at various angles of attack and sideslip. 

Major objectives of this study are: 1] to review some 

of the recent developments in dynamic distortion prediction 
with the Melick method as a foundation; 2] to demonstrate 
the utility of a new tool for improving the accuracy of peak 
instantaneous distortion contour maps; and 3} to evaluate 
present and recent developments in Melick dynamic distortion 
analysis . 
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1. REVIEW OF BASIC CONCEPTS 
AND DEVELOPMENTS 

A. The Melick Vortex Model 

The Melick convecting vortex model is a tool used to 
statistically determine the most probable peak instantaneous 
distortion level, given the steady-state distortion and the 
root mean square (rms) total pressure Fluctuation level at 
the engine face. It is formulated around the observation 
that the total pressure fluctuations exhibit random charac- 
teristics, with a near-normal (Beta/Gaussian) distribution 
(fig. 2]. From Bernoulli’s flow relationships, it is easily 
seen that these total pressure fluctuations can be expressed 
in terms of perturbations in the steady-state flow velocity* 
These velocity perturbat ions can in turn be modelled by 
time-variant vorticity (fig. 5]. Thus the Melick method 
envisions the total pressure fluctuations as being totally 
attributed to a series of random vortices (random in size ? 
strength, location, and orientation) convecting through the 
inlet duct (ref. 2). 

According to the Melick model, as a vortex passes 
through the inlet duct, it would create a fluctuation In 
the steady-state pressure level at all locations in the 
measurement plane, that is, the engine face. This pressure 
fluctuation would give rise to an instantaneous distortion 
level, computed from any of a variety of distortion factors 
(table 1). Given the properties of an arbitrary vortex, the 
resulting velocity perturbat ions can be determined from 
simple flow relationships (fig. 5). The pressure fluctuation 
can again be determined from the velocity perturbation, 
resulting in an instantaneous distortion level. 
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It is shown in reference 2 that the statistical proper- 
ties of the convecting vortices of the Melick model are 
directly related to the statistical properties of the pres- 
sure fluctuations. Specifically, the mean vortex size can 
be determined from the root mean square total pressure fluc- 
tuation level. This is accomplished by computing the rms 
fluctuation level resulting from an assumed vortex size and 
strength, and then comparing the measured rms level. The 
vortex size is then adjusted until the analytical and ex- 
perimental rms levels match. 

Once evaluated, the mean vortex properties are then 
used to compute the mean instantaneous distortion level, 
which leads to the determination of the most probable maxi- 
mum instantaneous distortion level. The mean instantaneous 
distortion is found analytically from the steady-state dis- 
tortion level and the rms total pressure fluctuation level, 
along with the mean vortex size (ref. 2). The peak instan- 
taneous value is then statistically extrapolated given the 
mean instantaneous value, the rms level, and certain statis- 
tical parameters (ref. 2). The maximum instantaneous distor- 
tion level can be computed for a variety of confidence 
levels, though the TT most probable” (a 50% confidence level) 
is used in most analyses (fig. 6). 

The newly computed maximum instantaneous distortion 
level is then used to produce the peak dynamic distortion 
contour map. First, the mean vortex is modified to accomo- 
date the peak dynamic distortion level. This is done by 
increasing the strength of the vortex until it produces 
an rms fluctuation level, and consequently a distortion 
level, which matches the maximum instantaneous distortion 
level. When this new vortex strength has been established, 
the resulting pressure disturbances are computed for each 
of the probe locations, and added to the steady-state pres- 
sures . The maximum instantaneous pressure array is then 
used to generate the peak dynamic distortion map (ref, 2 ). 
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B. The Minimum Dynamic Measurement: Approach 

One of The benefits of The Melick approach To dynamic 
disTorTion prediction is if T s low cosT relaTive To oTher 
meThods. TradiTionally , The Melick method requires steady- 
staTe ToTal pressure measuremenTs , along wiTh rms total 
pressure f lucTuaTion measurements aT forTy probe locations 
across The engine face. In The derivaTion of The mean in- 
stanTaneous disTorTion level, The mean value (face-average] 
rms level is used. The acTual number of high-response 
dynamic probes is noT imporTanT - jusT The mean rms value 
is of inTeresT. In principle, Therefore, The Melick method 
requires only one dynamic rms ToTal pressure measurement, 
provided an average value is indicaTed. In The inTeresT of 
furTher reducing insTrumenT aT ion cosT, and inleT blockage 
during a TesT run, if is desireable To minimise The number 
of dynamic probes used while refaining The accuracy of The 
resulTs. Proper placemenf of a minimum number of dynamic 
probes is necessary in order To obfain an accuraTe repre- 
senfafion of The average rms level, and The resulting peak 
dynamic disTorTion predicfion. 

Chen (ref. 2) provided a criTerion for The selecfion of 
dynamic probe locafions which yield reasonable accuracies 
in mean rms level deTermination. IT was observed That There 
exisTs an inverse relafionship befween The rms ToTal pres- 
sure f lucTuaTion level and The magnitude of The ToTal pres- 
sure. In oTher words, high rms pressure fluctuations Tend 
To occur in regions of low ToTal pressure , while low rms 
levels occur in high pressure regions. Furthermore if was 
noTed fhaf average rms levels Tend To occur near regions of 
average pressure. This implies fhaf dynamic probes placed 
near The steady-stafe mean pressure line would give rms 
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■total pressure fluctuation levels nearly equal to the face- 
average value. 

Since it is preferred to remain on the conservative 
side in dynamic distortion prediction, that is it would be 
more desireable to overpredict rather than underpredict the 
true peak dynamic distortion in any simplifications, it is 
suggested that the preferred dynamic probe location should 
be at or outboard of the mean pressure line (ref. 2). This 
will allow in most cases an rms level slightly higher than 
the average value obtained in a 40-probe analysis. In any 
case, dynamic probe locations selected should avoid regions 
of very high and very low steady-state pressures. 

The accuracy of this criterion is shown herein and in 
reference 2. It was shown that using 2 probes selected 
according to the "conservative side” criterion yielded 
distortion factor errors generally within 5% of the 40-probe 
prediction. Naturally, if dynamic probes were selected such 
that the average rms value were exactly equal to the 40- 
probe average, there would be no error. Conversely, the 
selection of improper probes can lead to very large errors. 
Consequently, the careful selection of locations for the 
placement of dynamic probes is extremely important for the 
accuracy of the results. 

C. The Turbulence Modelling Approach 

Because of the sensitivity of the predicted peak dis- 
tortion level to the indicated mean rms level, which in turn 
is sensitive to the location of the probes relative to the 
mean total pressure line, it is desireable to develop an 
approach which includes the benefits of both the full [ 40 - 
probe ] dynamic data method and the minimum dynamic data 
approach. In response to this need, Chen (ref. 3} developed 
a turbulence modelling approach which produces an accurate 
prediction of the peak instantaneous distortion with no 



8 


requirement for dynamic rms “total pressure fluctuation data* 

In this turbulence modelling approach, the rms total 
pressure fluctuation levels are simulated from information 
derived from the steady-state total pressure measurements* 
First, the axial velocity distribution (relative flow velo- 
city at each steady-state probe location] is calculated from 
the steady-state measurements. A set of turbulence modelling 
equations is then employed to compute the turbulent kinetic 
energy distribution, and the turbulent kinetic energy 
dissipation rate. These terms represent the turbulence 
levels required to generate the steady-state distortion. 

The rms total pressure fluctuation levels are then evaluated 
from the turbulent kinetic energy and the turbulent kinetic 
energy dissipation rate. These simulated rms levels are then 
used to compute the mean vortex properties, the mean Instan- 
taneous distortion level, and the peak dynamic distortion in 
the same manner as the original Melick model (ref. 3). 

The advantages of the turbulence modelling approach are 
obvious. There is no need for rms total pressure fluctuation 
levels to be measured - hence no high-response dynamic 
probes are needed. Instrumentat ion costs are reduced con- 
siderably from the fully instrumented 40-probe case. In 
addition, there is no need for concern over where to most 
effectively place a minimum number of dynamic probes. The 
turbulence modelling approach, when coupled with the Melick 
vortex model, is an efficient tool for determining the most 
probable peak instantaneous distortion level, given the 
steady-state measurements. 

The accuracy of the turbulence model is demonstrated in 
reference 3. It is shown that this approach is at least as 
accurate as the fully instrumented case in comparison to the 
□YNADEC results for a variety of inlet configurations and 
operating conditions (ref. 3). 
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2. A MODIFIED VORTEX MODEL 
A. Introduction 

Although the original Melick vortex modelling approach 
(including the modifications summarized in the previous 
section) is shown to be reasonably accurate in the predic- 
tion of peak instantaneous distortion levels, it is not as 
accurate as some methods in the generation of peak instan- 
taneous distortion maps (ref. 2, 5). It is therefore 

desireable to develop some modification to the Melick vor- 
tex model which can improve the accuracy of the peak instan- 
taneous map. 

It has been suggested (ref. 2) that the fault in the 
Melick peak instantaneous mapping method may lie in one of 
the vortex modelling assumptions. This modelling approach 
produces the peak instantaneous pressure distribution by 
superimposing a linear vortex along the mean shear line of 
the steady-state distortion pattern (fig. 3). The induced 
flow velocities produced by this vortex, whose properties 
are determined from the rms total pressure fluctuation 
levels, result in an amplification of the total pressure 
distribution. Both high and low pressure regions are en- 
hanced by this vortex so that the distortion level is magni- 
fied . 

In reality, the mean shear line of most steady-state 
distortion patterns is not a straight line, but is instead 
curved. In fact, often the mean line forms a distorted ring. 
This suggests that the core of the peak instantaneous vortex 
should not be a straight line, but should follow the curves 
of the mean shear line. This will be the basis of the pre- 
sent study. 



10 


B. Method of Approach 

In the present analysis, a new vortex modelling ap- 
proach designed to improve the accuracy of Melick peak 
instantaneous distortion maps is developed. In this new 
approach, the linear vortex model of the original Melick 
method (fig. 3] is replaced by a vortex which can have a 
nonlinear core (fig. 4). This is a justified modification 
because the mean shear line (the borderline between rela- 
tively high and low pressure regions) is generally non- 
linear ( f ig . 4 ) . 

There are three general methods in modelling a vortex 
with a nonlinear core. The first and most complex method 
would be to f or mu late a mathematical expression for a curve 
which fits the desired shape of the vortex core - that is 
the mean shear line. This expression could be in terms of 
2-dimensional cartesian or polar coordinates, derived from 
a least-squares (or other nonlinear) analysis, or perhaps 
from an infinite series expansion. This method has the 
potential of being extremely accurate as far as modelling 
the vortex is concerned, but would not be very efficient In 
terms of the computational effort. 

A second approach to modelling a nonlinear vortex core 
might be form a finite element model. The nonlinear vortex 
would be divided up into a series of linear vortex segments 
which would lie along the mean line. The number of segments 
used would depend on amount of curvature in the mean shear 
line and the desired resolution. This method, depending on 
the number of divisions selected, could be as accurate as 
the least squares/infinite series method, with considerably 
better computational efficiency. 

It is clear that these two methods have . the capability 



of achieving very high resolution in -the calculation of the 
vorticity effects, and in the generation of the peak instan- 
taneous distortion map. This high resolution capability is 
not necessarily useful, however. It should be kept in mind 
that the peak instantaneous map is generated by calculating 
the effect that the vortex has upon the total pressure 
readings obtained at the steady-state probe locations . Pres- 
sures at locations between probes are then interpolated from 
these new "readings”. Consequently any vortex action which 
occurs between probe locations is ignored, prior to the 
interpolation process. This limitation in the useable resol- 
ution of the vortex model is the basis of the simplifica- 
tions of the third modelling approach. 

The third approach to modelling a nonlinear vortex core 
is similar to the finite element model, but includes some 
important simplifying assumptions. First, the vortex is 
divided up into eight segments, each associated with one 
of the probe rakes [fig. 4). Each vortex segment is consi- 
dered the dominant contributer to the pressure disturbances 
occurring on the rake associated with that vortex segment. 

It is assumed that each vortex segment affects the pressure 
only on the rake associated with it. The position of the 
vortex segment relative to it*s associated rake is assumed 
to be at the probe nearest to the mean pressure line where 
it crosses the rake. The orientation of each vortex segment 
is assumed to be perpendicular to it’s associated rake, and 
coplaner with the measurement plane. Each of these simpli- 
fying assumptions are illustrated in figure 4, and are dis- 
cussed separately . 

The first assumption involves the division of the non- 
linear vortex into eight linear sub- vortices , or vortex 
segments. Traditionally, there are eight rakes mounted at 
45 degree intervals around the measurement plane (fig. 7), 
Since all probes on a rake are affected by the induced 
flow velocity caused by the local vorticity, it makes sense 
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■to divide the probes into rake-groups, and to determine the 
dominant vortex activity associated with that group. There- 
fore, in this study the nonlinear vortex system is divided 
into a set of linear vortex segments, with each segment 
acting as the dominant vortex activity for one of the rakes* 
Eight rakes each require one vortex segment, for a total of 
eight sub-vortices. Each vortex segment is considered by 
definition to affect only the probes on it T s respective 
rake, and induced vortex activity on adjacent rakes is con- 
sidered by definition negligible. 

The next assumption involves the definition of the lo- 
cation of each of the vortex segments. It is assumed that 
the vortex segment is placed directly over the probe loca- 
tion nearest to the mean shear line as it crosses over or 
passes near to the rake. In addition, it is assumed that 
the vortex segment is oriented perpendicular to the rake. 
These two simplifying assumptions are illustrated in figure 
3. It is suggested that these simplifications introduce 
only small errors into the analysis, while they allow con- 
siderable improvement in computational efficiency. In any 
case, the error produced by these simplifications will 
always be less than the error produced in the original 
Melick single linear vortex model. 

Finally, it is assumed that the vortex properties as 
derived in the Melick linear modelling approach are still 
valid in the segmented modelling approach. These properties 
include: 1) the mean vortex radius, and 2) the vortex 

strength. These terms were derived as a function of the 
rms total pressure fluctuation level, and the most probable 
peak instantaneous distortion level. 

Each of these assumptions and simplifications are made 
in the interest of providing a strai ghtf orward model and a 
simplified analysis. None of the assumptions are expected to 
introduce significant error into the analysis. The nature 
of the model and the analysis is intended to be preliminary. 
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in the interest of determining whether f urther 
warranted in this modelling approach. 

In the following section, the development 
matical formulations is presented based on the 
assumptions . 


research is 

of the mathe- 
simplifying 


C . Mathematical Formulations 

In the Melick approach to peak instantaneous distortion 
prediction, there are two distict sections: 1) the develop- 

ment of the most probable peak instantaneous distortion 
level; and 2) the generation of the peak instantaneous map. 
Since the present analysis is concerned primarily with the 
latter of these two sections, the first section will be 
presented only in summary form. Details on the derivation 
of the peak instantaneous distortion level may be found in 
reference 2. 

As described in section 1 . A . , the random total pressure 
fluctuations measured at the compressor face are attributed 
to the convection of a series of random vortices through the 
measurement plane. The pressure fluctuations are to be 
expressed in terms of velocity perturbations introduced by 
these vortices. The velocity profile of a one-dimensional 
steady and incompressible vortex is given as (fig. 5): 

v = V - e -^n(r/a) 2 - 1] Cl) 

T T a 
max 

where: V T is the tangential velocity at any radius r 
V 

T is the maximum vortex swirl inq velocity at r ; 

maX , . . I ^ 

a measure of vortex strength . 

r is the independant variable: radius 

a is the radius at the point of maximum swirling 
velocity - also called the vortex size 

e is the exponential 
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The total pressure Fluctuations produced by the vortex 
are superimposed onto the steady-state total pressure to 
Form a time variant instantaneous pressure: 

P T = P T + dP T (2) 

ss 

where: Pj is the instantaneous total pressure 

is the steady-state total pressure 
ss 

dP_j_ is the pressure Fluctuation produced by the 
vortex 

From the i ncompressible Bernoulli equation: 

P T = Pg + %P U 2 
SS 

p t = p s + %nu + v T ) 2 

where: P^ is the static pressure 

f is the Flow density 

U is the steady-state Flow velocity 

V T is the vortex- induced velocity 

2 

Let q = y 2 ?U be the steady-state dynamic pressure. Then 
substituting (3] into (2), we obtain: 

dP y = p y - p y = %Hu + m t ) z - y 2 ?u 2 

SS 

= y 2 f (u 2 + SUV ? + V 2 - u 2 ) 

2V T V T 2 2qV T 

= * c'-n 1 ’ + 'it 1 ] = nr 1 


(3a) 
( 3b j 


Second order terms have been neglected for V^. much less than 

U. 
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Substituting (1) into [4], we obtain: 


d p t =^ V T r e- ,S[tr/ » ) - 1] 

T Ua max 


cs: 


Equation (5) represents the total pressure fluctuation level 
produced by the convection of an arbitrary vortex through 
the inlet duct in terms the relative size and strength of 
the vortex, and the position of the probe relative to the 
vortex . 

In the Melick analysis, in order to determine the most 
probable peak instantaneous distortion level, the mean vor- 
tex size must be determined. This parameter is shown (Ref 2} 
to be a function of the rms total pressure fluctuations: 


dP„ 


dP 


rms ( f ] ^ 2 


T 


erf (7.38 £g) 


( 6 ) 


rms 


The quantity on the left hand side of equation (6] is the 
square of the ratio of the root mean square total pressure 
fluctuation level filtered at cut-off frequency f, to the 
unfiltered rms level. These quantities may be measured, 
or simulated as developed in reference 3. Using equation 
(B), the mean vortex size, a, can be solved for iteratively 
in terms of the filtered and unfiltered rms levels, the 
filter frequency, and the flow velocity, all of which are 
known quantities. 

The mean vortex size is then used to generate the mean 
instantaneous distortion level, k: (Reference 2] 

k = k + f ( a , dP ) (7) 

ss 5 T 


where: k = the steady-state distortion level determined 

gg ^ 

from the steady-state total pressure data and 
table 1 



Since dP^ is a function of the vortex properties, k can be 
be determined from them: 

k = f(P T ,q,U,S,V T ) ( 8 ) 

ss max 


where: , q, and U are measured quantities 

ss 

a is computed from equation [6] 


V^. is then determined, using equation (5): 

max 


UadPT VaCC r/a) - 1] 


T 2qr 

max 


( 9 ) 


For r = a, equation (9) is simplified to; 


max 


% - dP 

*3 1 rms 


( 10 ) 


Equation [8] then becomes: 

k = f(P T ,q,U,a,dP T ] [ 11 } 

ss rms 

Since q and U are constants, equations [7} and [11} are seen 
to be identical. 

The mast probable peak instantaneous distortion level 
is then statistically extrapolated as a function of the mean 
instantaneous distortion level, the rms total pressure fluc- 
tuation level, and a set of statistical parameters [ref. 2}: 

^max = ^Ck,dP_j_ , statistical parameters} [12J 

rms 


The vortex strength is then adjusted to match the change from 
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■the mean instantaneous distortion level to the maximum in- 
stantaneous distortion level: 


v t = V T + fCk - k) C i 3 ) 

max/pk max max 

Once the most probable peak instantaneous distortion 
level and the value of VTrnax has been determined, the 
value for the peak instantaneous total pressure can be de- 
termined : 


P 





[14) 


where dP is obtained from equations (5) and (13): 
pk 


dP 


T 


pk 


2q V T r e - 1 /aCCr/5) 2 

Ua max/pk 


1 ] 


C 153 


The only variable in equation (15) above is the radius, r. 

Using the simplifications and assumptions given in section 

2.B., the value of r can be determined on a rake-by-rake and 

probe-by-probe basis. Substituting (15) into (14) produces 

the following relationship (let V = V^_ ): 

max/ pk 


P (k, p) 
pk 


P T (k,p) + r~Ck,p) 


e -^[(r(k,p}/a) 2 


1 ] 


(IS) 


The subscripts k and p refer to the rake and probe number, 
respectively. The value of r(k,p) is defined as the distance 
between the probe with coordinates (k,p) and the core of the 
vortex associated with rake k: 


r(k,p] = r C k , p } - r v (k) 


C 17} 
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where: r^(k,p) is the radial location of probe (k,p) 

r Ck) is the radial location of the core of the 
v 

vortex associated with rake (k] 

The final step in generating the peak instantaneous map 

is to interpolate values for P at each of the discrete 

pk 

points between the probes. This is done in the same manner 
as with the steady-state map (ref, 6), 



3. RESULTS AND DISCUSSIONS 


Irrt *the Following sections, numerical and graphical 
predictions From the analytical methods described in the 
present work are provided with three inlet data sets. Data 
comparisons with the DYNADEC results are also provided with 
each oF the inlet conFi gurat ions . The three inlet data sets 
consist oF provisional experimental results From subsonic , 
transonic, and supersonic inlet conFi gurati ons under various 
Flight conditions. Inlet conFigurations and measured results 
oF the data sets are provided in Figures 8 through 10* Data 
comparisons oF predicted and measured peak instantaneous 
distortion levels, and graphical comparisons oF predicted 
and measured peak instantaneous distortion maps are also 
provided . 


A. Subsonic Inlet 

ConFiguration oF a Full-scale short S-shaped subsonic 
inlet duct is shown in Figure 8. The engine centerline is 
tilted approximately six degrees From the horizontal as 
shown in the Figure. The Freestream Mach number was given as 
subsonic. Six test cases were available For data comparison. 
These data were provided by the Air Force [AFWAL], Wright- 
Patterson AFB , Ohio. 

Comparison oF Melick predicted peak instantaneous dis- 
tortion levels is given in Table 2 and Figures 11 and 1 2 , 
Figure 13 shows mapping comparisons For the steady-state, 
DYNADEC measured peak, Melick predicted peak, and the Modi- 
Fied Vortex predicted peak instantaneous distortion patterns. 
As described in ReFerence 3, reasonably good accuracy of the 
distortion level prediction analyses is indicated. In certain 
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cases, -the peak instantaneous distortion level is underpre- 
dicted by the Melick approach. This is attributed to the 
fact that the Melick approach cannot accurately predict the 
peak distortion level far inlet flows with separated boun- 
dary layers. Unfortunately, the subsonic data set contains 
separated boundary layers at the engine face (ref. 3) Fur- 
ther study will be required to improve the Melick predictive 
accuracy in separated flow cases. 

The Modified ( segmented-nonl inear ) Vortex technique 
compares favorably with the Melick modelling approach in the 
peak distortion map generation, in certain cases. In cases 
where the Melick linear vortex model produces an accurate 
prediction of the peak instantaneous map, the modified 
approach generally overpredicts the distortion pattern 
slightly. In cases where the Melick approach yields poor 
results in the peak distortion map, the modified approach 
tends to improve the map considerably (fig. 13). 

B. Transonic Inlet 

Configuration of a 15% subscale long S-shaped transonic 
inlet duct is shown in Figure 9. Six test cases with a tran- 
sonic freestream Mach number were available for comparison. 
These data were provided by the Air Force (AFWAL) Wright- 
Patterson AFB, Ohio. 

Comparisons of predicted and measured peak instantane- 
ous distortion levels are given in Table 3 and Figures 14 
and 15. Figure 15 shows mapping comparisons for steady- 
state, DYNADEC measured peak, Melick predicted peak, and 
the Modified Vortex predicted peak instantaneous distortion 
patterns. Good accuracy in predicting peak dynamic distor- 
tion levels is indicated for these test cases. The Melick 
method slightly overpredicts the peak distortion level, 
which is the desired affect. 
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The Modified Vortex approach again compares favorably 
with the Melick linear vortex approach. In examples where the 
original Melick approach yields poor predictions of the peak 
instantaneous distortion pattern, the modified approach pro- 
duces superior results (fig. 16). In cases where the Melick 
approach produces good peak distortion maps, the modified 
approach produces comparable results. 

A notable exception to this can be seen in the first 
case (case number 464.12). The steady-state map shows a 
symmetrical pattern, while the DYNADEC predicted pattern 
is not symmetrical (fig. 16). The Melick approach produces 
a fairly symmetrical pattern as expected, while the modified 
approach produces a pattern almost identical to the steady- 
state pattern, except for enhanced pressure magnitudes, also 
as expected from the modelling criteria. This is due to the 
fact that this particular test case represents an extremely 
high angle of attack, where asymmetrical vortex shedding is 
evidently taking place. This inlet "pumping" has the effect 
of alternating high and low pressure levels on either side 
of the inlet duct instantaneously, while providing apparently 
symmetrical patterns in the steady-state. Asymmetrical or 
alternating vortex shedding is the same phenomenon which is 
associated with wing "rocking" in highly swept delta wings 
at high angles of attack. 

C. Supersonic Inlet 

Configurations for four 25% scale supersonic inlet 
ducts are shown in Figure 10. These inlet models include 
data for a variety of supersonic freestream Mach numbers, 
and angles of attack and yaw. There are thirteen test cases 
availabe for comparison. These data were also provided by 
the Air Force ( AFWAL ) Wright-Patterson AFB , Ohio. The four 
inlet conf i gurat 1 ons , test conditions, and some measured 
results are given in figure 10. 
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Comparison of predicted and measured peak instantaneous 
distortion levels is given in table 4 and figure 17. Figure 
18 shows comparisons of steady-state, DYNADEC measured peak, 
Melick predicted peak, and modified Melick predicted peak 
instantaneous distortion contour maps. Many of the peak 
instantaneous distortion levels are underpredicted slightly, 
primarily because these cases show separated boundary layers. 
It is recalled that the Melick approach tends to underpredict 
peak distort ion levels when separated boundary layers occur 
in the inlet duct. It is noted that in many cases the maps 
generated by the Melick approach appear to have no apparent 
pattern as far as relatively high and low pressure regions. 
The steady-state and peak instantaneous maps in these cases 
exhibit quasi-random character ist i cs , indicating severe 
turbulence levels and flow separation. The Melick prediction 
technique generally requires a reasonably well-defined mean 
shear line in order to effectively apply the vortex model. 

The predicted distortion levels and distortion maps are seen 
to be fair to good in these cases. In certain cases, when the 
linear Melick vortex model fails to provide a good prediction 
of the peak distortion map and the measured peak map resem- 
bles the steady-state pattern, the segmented vortex model 
provides a map superior to the one generated by the linear 
vortex model (fig. 13). 

D. General Results and Comments 
1 . The Melick Approach 

The Melick approach to predicting the peak instantaneous 
distortion levels can be evaluated by examining Tables 2, 3 

and 4, comparing peak distortion values as measured from the 
DYNADEC system, and as predicted by the Melick approach, 
Figures 11, 14 and 17 show these comparisons in graphical 

form . 
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It is seen from Figure 11 that: the subsonic peak dis- 
tortion factors are underpredicted in four out of six cases , 
and has a percent error of greater than plus or minus twenty 
percent in three out of six cases. At first glance this may 
be disturbing, but it is recalled that the subsonic inlet 
data set indicates a separated boundary layer, as described 
in Reference 3. Since the Melick approach assumes an attached 
boundary layer, the results are understandable. Nevertheless 
it can be said that the Melick approach provided a good 
ball-park figure in distortion level prediction. 

Figure 17 shows peak distortion level predictions with- 
in 20 percent in 4 out of 6 cases, with an underpredi ct i on 
of the measured peak distortion in only one case. These 
results can be considered very good. Near-perfect predictions 
are seen in two cases, which is encouraging. These transonic 
inlet data sets show a mean percent error of approximately 
ten percent for all six test cases. 

The supersonic test cases show peak distortion predic- 
tions within 20 percent in 11 out of 13 cases. However, it 
is also noted from Figure 17 that the peak distortion is 
underpredicted in almost all cases. As seen from Figure 18, 
the supersonic inlet cases in many cases represent highly 
turbulent separated flow conditions, for which the Melick 
technique is known to tend to underpredict. The overall re- 
sults for the supersonic test cases can be said to be fairly 
good, and very consistent. 

The overall accuracy of the Melick approach can be 
judged with Figure 19. The overall percent error in the 
predicted peak instantaneous distortion level with non- 
separated flow shows a mean value of +19.5%, while the 
separated flow inlets show a mean percent error of -3.3%, 
These results are considered good for preliminary engineering 


purposes . 
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2. The Minimum Dynamic Measurement Approach 

In addition to the results indicated in Reference 2 , 
which shows very good results in the minimum dynamic measure- 
ment approach with respect to the full ( 40-probe } approach , 
Tables 2 and 3, and Figures 12 and 15 indicate excellent 
correlations between the two approaches. In addition, as 
indicated in Reference 3, the turbulence modelling approach 
shows excellent predictions of peak instantaneous distortion 
levels, with no dynamic measurements. Predictions well with- 
in 20 percent of the 40-probe predicted values are indicated 
for most of the test cases. These results appear to validate 
these two low-cost approaches. 

3. The Segmented Vortex Approach 

The segmented vortex approach can be judged in terms 
of its performance with respect to theoretical expectations* 
The segmented vortex approach is, again, a simplified model 
of the nonlinear mean-line following approach. This model 
will always produce a predicted peak instantaneous distor- 
tion pattern similar to the steady-state pattern, with the 
pressure levels amplified somewhat. This phenomenon can be 
easily seen in the distortion map comparisons (figs. 13 , 16 

and 183. In cases where the peak instantaneous map Is not 
similar to the steady-state map, this approach will produce 
a poor prediction, while cases where the steady-state and 
peak instantaneous maps have similar patterns as measured 
by DYNAOEC, the approach will produce a map generally 
super ior to the linear Melick vortex approach . In some cases 
it is seen that the segmented vortex approach is far superior 
to the linear, while in most cases, the improvement is only 
marginal (fig. 13, 16, 18}. It is possible that the accuracy 

of this approach can be improved by removing simpl if Icat i ons , 
though the overall pattern would not change significantly. 
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CONCLUSIONS 


A simplified nonlinear vortex model has been developed 
in order to improve the quality of predicted peak instan- 
taneous distortion maps. The nonlinear mean-line following 
vortex model is simplified by dividing the vortex into 
linear vortex segments, one for each rake of probes, orien- 
ted perpendicular to each probe, and each having the charac- 
teristics of the mean vortex developed in the original Melick 
approach . 

A review and evaluation of recent developments in the 
Melick peak instantaneous distortion level prediction tech- 
nique has been included. A simplified description of the 
Melick method has also been provided, with references to 
more detailed reports. Predictions using limited, minimum 
and no dynamic data have been compared to OYNADEC measure- 
ments with favorable results for three inlet data sets. 
Minimum and no dynamic data approaches have also been com- 
pared to full ( 40 -probe) predictions with excellent results. 

The Melick approach, along with recent improvements, 
is shown to be an efficient and accurate design tool for 
predicting peak instantaneous disortion levels in prelimi- 
nary analyses. It is noted, however, that the approach does 
not work well with highly turbulent separated flow inlet 
conditions, due to limitations in the modelling approach. 
Further research will be required to develop improvements 
in separated inlet flow predictions. 

The segmented vortex approach is a useful method of 
improving peak instantaneous distortion maps, provided the 
Melick peak distortion level has been accurately predicted, 
and provided the actual peak distortion map does resemble 
the steady-state map pattern. Further improvements in this 
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modelling approach are possible by removing simplifications 
and assumptions, as long as these two conditions are met. 

At this time there exists no modification to the Melick 
vortex model which can accurately predict the peak instan- 
taneous distortion pattern when the measured peak pattern 
is significantly different from the steady-state pattern. 
Further study will be required to understand and predict 
this particular problem. 
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Table 1. Distortion Factor Definitions (Ref. 2) 


Factor 

Equation 

Supplemental equations 

Definitions 

I DC 

max 

IDC ^maxl^ClDC^IDC^) , 
max v 1 2 

»5(IDC a +IDC 5 )^ 

(p^) 4 ~(p^ . ), 

IDC - LJ t 1 

3 p t 

(P t )j“ average total pressure for ring j 

(P t m i n ^ 4 = minimum total pressure reading 
’ J in ring j 

P t ~ average total pressure at engine face 

IDR 

max 

IDR = max(IDR, ,IDR C ) 
max 1 5 

P t " 

IDR - LJ. 

J Pt 

K D2 

NR 

2 e.(Ap/p )(OD/D.) 

„ _ J=1 J C J 

H02 NR 

Z (OD/D.) 
j = l J 

4p ‘ . <Vr<Pt..i.> 1 4100 

p t <p t >j 

( Pt ) j» (p t,min ) j = see above 
0.** circumferential extent of largest 
continuous total pressure depression 
below (p^)^, degrees 

Dj- diameter of ring j ; NR ■» number of 

OD ® outer duct diameter r i n 8 

K e 

RR 

Z, (Al) (1/D.) 

Kn= J— J i— 

6 _ NR 

(q/p )Z (1/D ) 
j = l J 

<Vj- < a ; +b i >j 

<■.>,- b® |f -(v) 
‘Vj- sfe If ,ln<9 i0 1 

4p t1 ‘Pc’ 1 ^ p t .base' 1 

P t P t 

P t »(p t ) see above 

q - average dynamic pressure at engine 
face 

M * number of rakes 

(Pi.), 35 individual total pressure; rake i 
J , ring j 

0 « angular position of p t ^ 

— *-= r = base radial profile for 

P t ring j; set = 1 for all j 

b = radial distortion weighting factor=1.0 

k rad 

NR A Pt p 

Krad p t q ^ 

K A2 

K A2- V bK RAD 

5spr 

ASPR ASPR 


p^= see above 

p = minimum total pressure at engine 

t- * lliX 11 

face 

k = compressor reduced frequency 

Dist 

D1St Warnin' /p .) " 

ID 

ID = k (IDC)b + k (IDR) 

c r 


k^~ circumferential distortion sensitivity 

factor 

k = radial distortion sensitivity factor 
b « circumferential distortion weighting 
factor 










Figure 2. Illustration of a Typical Inlet Test Model 
and Peak Distortion Factor Measurement 
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Figure 8. Illustration of a Subsonic Inlet Test Model and sane lest 
Results (unpublished data from Air Force Flight Dynamics 
Laboratory, Wright-Patterson Air Force Base, Dayton, Ohio) 
[Ref. 3] 




























Figure 9 . Illustration of a Transonic Inlet Test Model and seme Test 
Results (unpublished data from Air Force Flight Dynamics 
Laboratory, Wright-Patterson Air Force Base, Dayton, Ohio) 
[Ref. 3] 
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ORIGINAL PAGE 18 
OF POOR QUALITY 


(a) Configurations of four .25 scale Tailor-Mate Model: 



Figure 10. Illustration of four Supersonic Inlet Test Model and some 
Test Results (ref. 3 } 
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Table 2. Subsonic Inlet Distortion 
Factor Comparison 
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Figure 11. Comparisons of the Predicted and Measured rms Level and 

Peak Distortion Factor for the Subsonic Inlet Model shown 
in Figure 8. (unpublished data) 
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IDC - 40 probe prediction 

max 


Figure 12. Comparisons of peak distortion factors predicted by the present 
analyses and the Melick method based on the total pressure rms 
measurements for the subsonic inlet data set 
(unpublished data) 
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Figure 13b. Subsonic Inlet Map Comparison [54.30] 
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Figure 13c. 


Subsonic Inlet Map Comparisons [81.40] 
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Figure 13d. Subsonic Inlet Map Comparisons [111.30 
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Figure 13e. Subsonic Inlet Map Comparisons (112.30] 
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(measured) K^-, 

m 


Figure 14. Comparisons of the Predicted and Measured rms Level and 

Peak Distortion Factor for the Transonic Inlet Model shown 
in Figure 9 . (unpublished data) 
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Figure 15. Comparisons of peak distortion factors predicted by the present 
analysis and the Melick method based on the total pressure rms 
measurements for the transonic inlet data set 
(unpublished data) 
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Figure 16b. Transonic Inlet Map Comparison (465.11) 


53 


percent from average pressure 



k-theta 

K02 

(IDC)-MAX 

(IDR)-^AX 

kra 

KA2 

OSPrt 

ID 


0.27191 

7ad.a 

0.04194 

0.05092 

0.17237 

0.40705 

0.03074 

0.73617 




0.440^2 

1157.7 

0 . 0-3070 

0. to925 
0.54299 
0.04502 
0 . 842o6 


steady- state 


□YNADEC peak 


percent from average pressure 



K-THETA 

K02 

( IOC 0 -MAX 
(I0R5-MAX 
KRA 
KA2 
DSPR 
10 


0.37410 

1230.5 

0.06921 

0.05092 

0.17237 

0,50924 

0.04214 

0.85062 



t IUC3-MAX 
ll jtf ; -HAX . 
KPA 
K 42 
L-5PP 
— 114.., 


0. U5«&9 
0.06169 
4,21/20 
0 . 5u660 
0.0^6 16 


Melick model peak 


segmented vortex peak 


Figure 16c. Transonic Inlet Map Comparison (473.12) 
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Figure 16d, Transonic Inlet Map Comparison [485.10) 
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Figure 18b, Supersonic Inlet Map Comparisons (189/3) 
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Figure 18c. Supersonic Inlet Map Comparison (216/3) 
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Figure 18e. Supersonic Inlet Map Comparison [246/3} 
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Supersonic Inlet Map Comparison [247/23 
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Supersonic Inlet Map Comparison ( 640 /S) 
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Figure 18h 


Supersonic Inlet Map Comparison (643/3] 
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Figure 18 j. Supersonic Inlet Map Comparisons ( 1334/2 
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Figure 18k. Supersonic Inlet Map Comparison (433/3) 
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Figure 181 


Supersonic Inlet Map Comparison (473/3) 
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Supersonic Inlet Map Comparison (1554/4) 




Figure 19. Accuracy of the Present Method in Predicting rms Levels 
and Dynamic Peak Distortion Factors 






PART II 


COMPUTER PROGRAM DOCUMENTATION AND USER'S MANUAL 

ESTIMATING MAXIMUM INSTANTANEOUS INLET FLOW 
DISTORTION FORM STEADY-STATE TOTAL PRESSURE 
MEASUREMENTS WITH FULL, LIMITED, OR NO DYNAMIC 

DATA 



ABSTRACT 


A computer program for statistically predicting peak 
instantaneous dynamic distortion, given steady-state distor- 
tion data and dynamic root mean square pressure fluctuation 
levels in gas turbine inlets, is presented. The statistical 
approach utilizes a physical flow model which character izes 
inlet flow distortion as due to random vorticity convecting 
through the inlet duct. Characteristics of a mean vortex are 
statistically determined to match steady-state distortion 
data and contour map, as measured by steady-state total 
pressure probes. The mean vortex character ist ics are then 
intensified according to the mean rms fluctuation level as 
measured by full or limited high response pressure trans- 
ducer instrumentation, or as simulated by turbulence model- 
ling, to produce the most probable peak instantaneous dis- 
tortion level. The computer program utilizes this approach 
to solve for the dynamic distortion and print the results, 
including contour maps. 
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FOREWORD 


This Report is designed to be a User r s Manual and Doc- 
umentation Guide For The improved Melick (ref . 1-3} dynamic 

distortion computer program developed at the University of 
Kansas. This program character izes the random vortices used 
to describe the unsteady, turbulent flow in jet engine inlets, 
and statistically calculates the most probable peak instan- 
taneous (dynamic) distortion level for a particular inlet 
operating condition. Steady -state distortion levels are com- 
puted for eight common distortion factors given the time- 
averaged steady-state probe pressure array , and the root 
mean square pressure fluctuation levels are used to project 
the maximum peak instantaneous distortion for the given con- 
ditions . 

Details of the derivation and development of the random 
vortex modelling approach are not included in this User f s 
Guide as the Guide is orientated more towards application 
than theory. The References, however, provide exhaustive 
detailing of the general Melick approach, especially Refer- 
ences 1 through 3. Reference 4 provides an extensive list of 
other sources which relate to distortion prediction. Finally 
References 5, 6, and 8 show details of some specific devel- 

opments in distortion research at the University of Kansas. 

Major segments in this Guide include descriptions of 
the main program and subprograms as well as input and output 
data. Sample problems are included for i 1 lustrat ion . A list- 
ing is provided in an appendix. The fully documented program 
requires memory capacity for 70,000 characters in 2300 lines. 
Hardware requirements include, in addition to the mainframe 
computer, an on-line printer for high speed output. 
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of dynamic data input or turbu- 
lence modelling subroutines 

30 

P 

PI 

Local instantaneous total pres- 
sure recovery 

28 

PAVG 

PAVE 

Average pressure recovery at 
compressor face measurement plane 

37 

PFIX 


Subroutine - transforms input 
data into mapping parameters 

13 

PFX 

PFXL 

Subroutine - computes vortex 
flux rate 

13 

PRNT 


Subroutine - prints pressure 
arrays and other output 

13 

PROBE 

NPR, NT 

Numerical designation of probe 

31 

PS 


Local steady-state total pres- 
sure recovery 

28 

PSAVG 

SMAVG 

Average steady-state total pres- 
sure recovery at engine face 

37 

PSD 


Power spectral density function 

- 

PSI 


Aicraft/ inlet yaw angle, degrees 

29 

P SPEC 

K, KD, OK 

General flow distortion factor 
designator 

26 

PT 

PTS 

Local total pressure recovery 
at compressor face 

37 

PTAVG 

PTAV , TMAVG 

Face average total pressure 
recovery 

37 
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OTHER NAMES DESCRIPTION 


PAGE 


PTMAX 

PTMX, TMMAX 

Maximum indicated local total 
pressure recovery 

37 

PTMIN 

PTMNj TMMIN 

Minimum indicated local total 
pressure recovery 

37 

QAVG 

QAV 

Face average dynamic pressure 

38 

QPT2 

Q/PT2 

Ratio of dynamic-to-total pres- 
sure at engine face 

42 

RAOLOC 


Radial position of rings, deg. 

25 

RATK 


Subroutine determines the effect 
of filter frequencies on dis- 
tortion factor values 

14 

R I 


Inlet centerbody hub radius in 
inches (typical units) 

41 

RKMN 


Average total pressure recovery 
along rake with minimum average 
pressure 

39 

RKMX 


Average total pressure recovery 
along rake with maximum average 
pressure 

40 

RKP 

KR 

Radial weighting factor for 
solution of ID distortion 

30 

RMS 


Root mean square 

- 

ROUT 


Percent difference between the 
indicated local probe , rake, or 
ring pressure and the face 
average pressure 

88 

RRT 

R/RT 

Radial location of vortex core 
(nondimensional) 

49 

RS 


Ratio of filtered-ta-unfiltered 
rms total pressure fluctuations 

32 

RSIGMA 


Subroutine - computes rms value 
of distortion factors 

14 
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SYMBOL 

RT 

SEARCH 

SG 

SGDK 

SGP 

SG/PT2 

SMAVG 

SPTRC 

SUMMER 

SYMBLE 

T 

TOPI 

TDP2 

TDP3 

TE 


OTHER NAMES DESCRIPTION 


Outer inlet: radius (tip) at the 
compressor face, in 

Subroutine - controls solution 
of peak distortion factors 

SIG Unfiltered rms total pressure 

fluctuations 

Filtered rms distortion level 

Unfiltered rms distortion level 

Ratio of unfiltered rms total 
pressure fluctuations to aver- 
age total pressure recovery at 
engine face 

PSAVG Face average static pressure 

recovery 

Total pressure recovery through 
supersonic inlet shock system 

Subroutine - evaluates error 
function in vortex core size 
iterative computations 

Subroutine - supplies symbols 
for distortion contour map 
generation 

TOP Time-on-point or data point 

dwell time - duration of rms 
pressure fluctuation data meas — 
urement , seconds 

Simple distortion parameter: 

( PTMAX-PTMI N ) /PTMAX 

Simple distortion parameter : 
(PTMAX-PTMIN) /PTAVG 

Simple distortion parameter: 
(PTAVG-PTMIN)/PTAVG 

E Turbulent kinetic energy dissi- 

pation rate - ft 2 /sec 3 
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41 

14 

32 

46 

46 

45 

37 

30 

14 

15 
27 

95 

95 

95 

45 
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SYMBOL 

OTHER NAMES 

DESCRIPTION 

PAGE 

THE 

THETA 

Angular location of the center 
of the arc DTH , in degrees 

40 

THMN 


Angular location of the rake 
with minimum rake-average total 
pressure recovery , degrees 

39 

THMX 


Angular location of the rake 
with maximum rake-average total 
pressure recovery , degrees 

39 

TITLE 1 


Identification block, also 
serves to execute end of file 

23 

TITLE2 


Identification title block #2 

28 

TMAVG 

PTAVG 

Average total pressure recovery 

37 

TMMAX 

PTMAX 

Maximum local total pressure 
recovery 

37 

TMMIN 

PTMIN 

Minimum local total pressure 
recovery 

37 

TRIDIA 


Subroutine - solves tridiagonal 
matrix equations 

15 

TURBUL 


Subroutine - manages turbulence 
modelling predictions and other 
computations as a subdriver 

16 

TV 

THETA 

Angular location of vortex core 
in degrees 

49 

UNSTDY 


Subroutine - manages dynamic 
data and peak distortion pre- 
diction as a subdriver 

15 

UU 

TUU 

Sum of squares of radial and 
circumf erential velocity gradi- 
ents 

45 

U 


Dimensionless velocity at each 
probe location 

76 

112 


Face average flow velocity , fps 

29 
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OTHER NAMES 


DESCRIPTION 


PAGE 


VBAR 



Average vortex strength, in^/sec 

49 

VBM 

VBMX 


Maximum vortex strength, in^/sec 

SO 

VBO 



Approximated vortex strength 

50 

VL 



Vortex core length limit 

49 



(Nornenclature Addendum) 
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ANRTU 



Function subprogram - computes 
vortex Flux rate 

7 

KD 

• 
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primary distortion Factor 

as 

NPG 



Data run i dent i Ficat ion code, or 
”part-point" 

31 
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I. INTRODUCTION 


Turbulence and other Flow nonuniformities in aircraft 
engine inlets have long been known to cause an unwanted Flow 
distortion phenomena at the compresson Face. These imperfec- 
tions in the ideally smooth inlet airFlow is Frequently due 
to the turning and shaping oF the Flow as it passes through 
the inlet duct. Generally the magnitude oF the distortion is 
a Function oF the angle oF attack and sideslip (yaw angle) 
of the aircraFt . The time averaged steady-state distortion 
level is relatively easy to determine experimentally by lo- 
cating an array oF total pressure probes ahead oF the com- 
pressor Face, and evaluating specific distortion parameters 
based on these steady-state measurements. Steady-state dis- 
tortion can be of suFFicient magnitude to disrupt the proper 
operation oF the engine by stalling the compressor. EFForts 
to develope high perFormance engine and inlet conFigurations 
has been hampered because oF the inherent sensitivity of 
highly loaded compressors to Flow distortion. 

It has also been Found that random Fluctuations In the 
distortion level, known as dynamic distortion, can have an 
even greater effect an engine stability as the steady-state 
distortion. It has been demonstrated that the dynamic dis- 
tort ion can cause the engine to surge even though the steady 
state distortion is well below the level at which the engine 
would be expected to stall. It becomes of particular impor- 
tance to be able to predict the dynamic distortion levels 
which could occur at any instant in time. 

One of the most common experimental methods of deter- 
mining the maximum i nstantaneous distortion is to use Fast- 
response (dynamic) probes to produce time histories of the 
total pressure Fluctuations at the compressor Face. These 
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instantaneous pressures are "then translated, as in the steady 
state case, into distortion parameters . These data are then 
screened by the Dynamic Data Editing and Computing System, 
DYNADEC , to determine the maximum instantaneous distortion 
during the test run. An estimation of the most probable peak 
distortion level is then available for the inlet designer . 

The DYNADEC approach to dynamic distortion prediction 
is generally quite accurate , but is extremely expensive in 
terms of test instrumentati on and computing time. For pre- 
liminary design purposes , it becomes difficult to justify 
the cost of a full DYNADEC test run. It is for this reason 
that methods of statistically predicting peak distortion 
levels have been developed. Further information on the 
DYNADEC and various statistical prediction methods can be 
obtained with the aid of References 4 and 7. 

Of the many statistical approaches for predicting peak 
dynamic distortion, the Melick random vortex model (Ref. 1 - 

3] is of particular interest because of it’s high efficiency 
in terms of data requirements and numerical analyses. The 
basis of the Melick approach is formulated around the obser- 
vation of the randomness of the total pressure fluctuations 
during a test run. It was hypothesized that the inlet flow 
could be modelled as having randomly distributed vortices of 
random strength, size, and orientation convecting with the 
steady-state flow, which itself is distorted by a large 
steady-state vortex. By applying fluid mechanics to convsc- 
ting vortices, a mathematical model of the inlet turbulence 
can be generated. The vortices are then translated into dy- 
namic distortion parameters using a statistical criterion. 

The distortion level, or the extent to which the flow 
is distorted, is generally defined in terms of distortion 
factors. These distortion factors are designed to indicate 
the distortion relative to some reference value, typically 
the level at which the engine could be expected to surae, 
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The maximum dynamic distort: ion prediction in the Me lick 
approach makes use of rms total pressure fluctuation levels 
to identify the main variables in the convecting vortex flow 
model [Ref. 5]. Filtered and unfiltered rms levels are re- 
quired so that any unwanted effects, such as engine speed, 
can be removed. The rms levels are somewhat easier to pro- 
cess than the instantaneous distortion computations done by 
□YNADEC , but the instrumentation requirements are much the 
same. It is seen that 1 nstrumentat i on costs can be reduced 
by using fewer dynamic probes . The Melick method allows a re- 
duction in probes since it actually uses the face-averaged 
fluctuation level in the analysis. In principle, the use of 
very few dynamic probes is feasible, as long as they produce 
the same face -average rms fluctuation level as the fully in- 
strumented case. 

There is some difficulty in choosing the locations for 
the placement of a limited number of dynamic probes, because 
it requires some knowledge of the solution before the test 
begins. Reference 5,- however, provides a simplified scheme 
for locating as few as two dynamic probes at the engine face 
while retaining sufficient accuracy in the dynamic distortion 
prediction. It is apparent that even further cost reductions 
could be achieved if the requirement for dynamic probes and 
the associated instrumentation could be eliminated entirely. 
Until recently, however, no methods have been available for 
reasonably accurate peak dynamic distortion prediction with- 
out dynamic data. 

Research at the University of Kansas has produced a 
new technique for estimating maximum instantaneous distortion 
based only on the steady-state total pressure measurements. 
Chen (Reference 6} has developed an approach to inlet turbu- 
lence modelling which analytically simulates the rms total 
pressure fluctuation levels using the predicted turbulent 
kinetic energy distribution at the compressor face. These 
simulated rms levels replace the rms level data which had to 
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The simulated rms fluctuation levels 


be measured previously, 
are then used to compute the variables of the random vortex 
model, from which the peak dynamic distortion parameters are 
derived, just as if the rms levels had been input as data. 

The purpose of this work is to present a computer pro- 
gram which statistically computes the most probable peak dy- 
namic distortion level, based on the methods of Chen and 
Melick. The program is designed to be highly adaptable in 
that the user may decide on the extent of the dynamic data 
to be input. There are three main alternatives available to 
the user. First the user may select a full set of dynamic 
data, a partial set can be considered (to a minimum of two 
dynamic probes], or the user can opt to input no dynamic da- 
ta. In the last case, the program automatically executes the 
computations related to the turbulence modelling and dynamic 
data simulation. This flexibility is designed to not only 
allow the user to select and contro 1 the quantity of dynamic 
data to be processed, but to also allow comparison of dif- 
ferent dynamic probe configurations in a single data run. 

In summary, the subject computer program solves for an 
estimation of the maximum instantaneous distortion, given 
the steady-state distortion data and rms total pressure fluc- 
tuation levels. The mathematical and theoretical derivations 
are well documented in the Melick references (Ref. 1 - 3 ) 
and the improvements by Chen are detailed in References 5 
and 6. Additional information on inlet flow distortion In 
general can be found with the aid of Reference 4. 

This Users Manual is designed to assist the user toward 
an understanding of the operational capabilities of the pro- 
gram . The three major sections of the Manual include a break- 
down of the program elements, an input and output data sec- 
tion, and a set of sample problems. A listing of the program 
is included at the end of the Manual. Suggestions for passi- 
ble future studies in improving the program or the analytical 
techniques are also included. 
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II. PROGRAM DESCRIPTION 


The subject computer program, the MAXIOYN peak dynamic 
distortion estimator, is written in FORTRAN IV and can be 
run as is or with minor modification on most FORTRAN compil- 
ers. MAXIDYN requires memory capacity for about 70,000 words 
in 2300 program lines, including comments. Deletion of the 
comment lines would reduce the memory space needed to about 
45,000 characters and 1600 lines. Appendix A of this work 
includes a listing of the program and subprograms. 

This program is designed to be flexible in nature, and 
can be used to run with a variety of inlet pressure probe 
conf igurations . Individual test cases can be analyzed sepa- 
rately, or groups of data sets can be run in sequence. A set 
of typical distortion factors are included in the program, 
though these can be modified by the User. Figure 7 gives the 
definitions of the distortion factors used in this program. 

The program may be used with or without dynamic rms 
total pressure fluctuation data, with a minimum of two probes 
in the case that dynamic data is included. 

A block diagram of calling sequences of subprograms is 
given in figure 1 . A description of each of the subprograms 
is given in this section. The subprogram descriptions are 
alphabetized, for convenience. An operational sequence of 
events is included to illustrate key events during a data 
run . 

Peripheral requirements are limited to a line printer. 
The program is suitable for use online while data is being 
collected, provided format requirements (sect. Ill) are met. 
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II. a. SUBPROGRAM DESCRIPTION 


1 . MAXIDYN main driver 

The main driver of The MAXIDYN program controls some 
of the data input, including the inlet probe configuration , 
and the steady -state pressure array. In addition, the main 
driver controls the subprograms which handle the remaining 
input data, distortion computations, and the output. Speci- 
fically, the main driver controls directly the following: 

a. Reading in of pressure probe ring and rake geometry. 

b. Reading in of data titles and identifying comments 
after checking for an End Of File command which 
stops program execution. 

c. Reading in of steady-state pressure data. 

d. Controlling the subroutines which control other da- 
ta input, check far errors, assign default values, 
and control distortion computations and output. 

2. Function ARNTU 

This function subprogram computes the vortex flux rate 
and it r s effect on the root mean square distortion level. 
ARNTU is controlled by subroutine MAXDP. 
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Subroutine CUBIC 


3 . 


This subroutine controls the cubic spline interpola- 
tions for subroutines TURBUL and INIVEL. These slope-based 
cubic spline interpo lat i ons are used to compute velocity 
gradients and turbulent kinetic energies at the fine grid 
points during turbulence modelling computations. 


4. Subroutine DISPAR 

This subroutine is used to calculate the eight distor- 
tion factors used by the program. These distortion factors 
are defined in figure 7 and can be modified by the user, 
DISPAR is controlled by subroutine DISTRT, a subdriver which 
controls most of the distortion computations. The actual 
formulas for the distortion factors are contained in DISPAR. 


5. Subroutine DISTRT 

Subroutine DISTRT is a subdriver which controls the 
computation of the distortion factors. Some of the duties of 
DISTRT includes the following: 

a. Calculation of simple distortion parameters ; for 
instance, the locations of the rake or ring with 
maximum and minimum average pressure . 

b. Calculation of average static pressure at the engine 
face, and the average Mach number. 

c. Control subroutines INTERP and DISPAR which continue 
the distortion factor computations. 
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6 . 


Subroutine EXTRME 


This subroutine manages the computation of extreme val- 
ues of the distortion factors. Called by subroutine MAXDP , 
EXTRME controls the solution of the most probable peak dis- 
tortion level for each of the distortion factors. The peak 
distortion factor is calculated by adding an incremental dis- 
tortion level to the steady state distortion. The incremental 
distortion level is computed via the SEARCH subroutine. tX i RME 
returns the peak distortion level to MAXDP after summing the 
steady state and incremental distortion values. 

7. Subroutine FINITE 

FINITE is a subroutine which is used to solve the fi- 
nite difference equations for subroutine TURBUL . These equa- 
tions are the turbulence modelling set formed by an implicit 
tridiagonal matrix scheme. The elements of the tridiagonal 
matrix equations, which consist of the turbulent kinetic en- 
ergies and the turbulent kinetic energy dissipation rates, are 
formed by FINITE and solved by subroutine TRIDIA. FINITE also 
computes the relative errors in the turbulent kinetic energy 
and the tubulent kinetic energy dissipation rate for each of 
the fine grid points at the compressor face. 

8. Subroutine FRF 

This subroutine evaluates the mean vortex core size by 
an iterative inverse solution scheme. FRF evaluates the vor- 
tex core size as a function of the f i Itered-to -unf i ltered 
root mean square total pressure fluctuation level. Subroutine 
SUMMER evaluates the error function of the vortex core size, 
and the solution is iterated until the error is small. 
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9 . 


Subroutine INITL 


This subroutine solves for the initial values of the 
turbulent kinetic energy and the turbulent kinetic energy 
dissipation rate. These initial values are used as a starting 
point in the iteration of the solution of these parameters. 
INITL is controlled by the TURBUL subdriver, which uses the 
turbulent kinetic energies to solve far the turbulent model 
in the synthesis of the rms pressure fluctuation levels. 


10. Subroutine INIVEL 

This subroutine calculates the circumferential and ra- 
dial velocity gradients at each of the grid points at the 
compressor face. INIVEL is called by TURBUL and uses subrou- 
tine CUBIC to carry out spline interpolations of the velocity 
gradients . 


11. Subroutine INTERP 

This is an i nterpo lati on subroutine which calculates 
the total pressure recovery at each of the discrete paints 
in the measurement plane. These points are used to generate 
the pressure contour map. INTERP uses linear interpolation 
to find the pressure at points between the pressure probe 
locations. Two linear interpolations are carried out: a ra- 

dial one and a circumferential one. The final value is taken 
to be the average of these i nterpo lat ions . INTERP is called 
by bath MAINLP and DI5TRT; when called by MAINLP, the Inter- 
polated values are used to generate the contour map, while 
DISTRT uses the interpo lati ons to compute the distortion fac- 
tors. A call to subroutine PRNT has been nulled - it had pro- 
vided a message when i nter po 1 at i o ns could not be performed. 
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12 . 


Subroutine LNPOUT 


This subroutine controls some of the output from the 
program . When called by UNSTDY , LNPOUT prints two of the ta- 
bles in the output: the Overall Flow Descriptors and the Flow 
Distortion Factors. The Overall Flow Descriptors table gives 
values for some of the simple distortion parameters , and the 
Flow Distortion Factors table gives values for the eight 
user-defined distortion f actors . LNPOUT prints these tables 
for both the steady-state and peak instantaneous case. In 
addition, LNPOUT prints the distortion contour maps for the 
steady-state and peak instantaneous case. UNSTDY controls 
LNPOUT by passing a control parameter; LNPOUT then selects 
the output to be printed. 

13. Subroutine MAINLP 

Subroutine MAINLP controls the calculations involved 
in the development of the pressure distortion contour map. 
MAINLP calls on INTERP to calculate the pressure at any of 
the discrete points at the compressor face, given the pres- 
sure at the probe locations. Subroutine SYMBLE then assigns 
a symbol for each of the discrete points, based on the pres- 
sure found by INTERP. MAINLP then passes the pressure and 
symbol information to the main driver, and ultimately to 
LNPOUT for printing of the distortion map. 

14. Subroutine MAXDP 

This subroutine is a subdriver which controls the com- 
putation of the peak distortion levels for the eight distor- 
tion factors. MAXDP computes the mean vortex size and the 
mean rms pressure fluctuation level, from which the peak 



instantaneous distortion is derived . MAXDP also controls "the 
calculations involved in producing the effects of other para- 
meters on the peak distortion, like the vortex flux rate and 
and engine filters, via subroutine RATK . Subroutine RSIGMA 
is called to compute the filtered and unfiltered rms level 
of the distortion f actors . Subroutine EXTRME then computes 
the peak instantaneous distortion statistically at 50%, 95%, 

and 99.7% confidence levels. MAXDP then prints the results 
in tables, namely the Distortion Factor Extreme Value table 
and others. The user-selected most probable peak instanta- 
neous distortion factor is also printed. It is this distor- 
tion factor that the peak distortion map is based in the 
iterative matching process. 


15. Subroutine NEWPSD 

This subroutine is a major subdriver which controls 
some of the input data and manages most of the computations 
involved in the prediction of the peak instantaneous distor- 
tion. NEWPSD controls the input of the dynamic data, includ- 
ing the filter frequencies, the datapoint dwell time (time 
on point], identification and program control parameters , 
and the rms pressure fluctuation data at each of the dynamic 
probe locations. NEWPSD also passes program execution to the 
turbulence modelling subroutine, TURBUL , if the user has se- 
lected the option of not entering dynamic data. Once the dy- 
namic data has been entered or synthesized, NEWPSD controls 
subroutines MAXDP and FRF which manage the computations in 
the peak instantaneous distortion prediction. NEWPSD also 
prints the dynamic data and the i dent i f i cat i on and control 
parameters . NEWPSD is controlled by the main driver and sub- 
routine UNSTDY 
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16. Subroutine PFIX 

This subroutine has two primary functions. The first 
part of the routine transforms the steady-state or dynamic 
distortion data from pressure recoveries to percent differ- 
ences from the average values. These percents are then used 
by the mapping routines for the plotting of the pressure 
distortion contour maps. PFIX also calculates the dynamic 
pressure and Mach number at each of the pressure probe loca- 
tions as a secondary function. The face-average Mach number 
is also computed by PFIX. PFIX is called by the main driver 
in the steady-state case, and UNSTDY in the peak dynamic 
estimation case. 


17. Subroutines PFX and PFXL 

These twin subroutines are used In the computation of 
the eddy [vortex) flux rate as a function of the distortion 
level. The difference between the two subroutines is in the 
computation of the vortex flux rate which depends on the 
magnitude of the ratio of steady-state to root mean square 
distortion: when this ratio is greater than 2.0, PFX is 

called, while PFXL is called when the ratio is less than 2 . 
The computational procedure for these two cases is somewhat 
different and an error would probably occur during computa- 
tions which involve logarithms and exponentials if the cases 
were not separated. 


18. Subroutine PRNT 

This subroutine controls the printing of steady-state 
and peak instantaneous pressure arrays, the printing of some 
of the titles and the listing of messeges in the output. 
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19. 


Subroutine RATK 


RATK is a subroutine which evaluates the effect of the 
engine filter frequency, FO, on the root mean square distor- 
tion level. The variation of rms distortion with engine fil- 
ter frequency is analytically determined. 


20. Subroutine RSIGMA 

This subroutine is called by MAXDP and computes the 
root mean square distortion level for the eight distortion 
factors. The routine is divided into separate groups for 
individual distortion factor evaluations. 


21 . Subroutine SEARCH 

Subroutine SEARCH controls the computation of a peak 
distortion parameter which is used by EXTRME to form an es- 
timation of the peak instantaneous distortion level. The 
ratio of the difference between the peak and steady-state 
to the rms distortion is solved for in an iterative search 
for the peak distortion level. 


22. Subroutine SUMMER 

This subroutine evaluates the error function in the 
iterative calculation of the mean vortex core size. SUMMER 
is called by FRF. 
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23, Subroutine SYMBLE 

This subroutine supplies the mapping symbols for the 
generation of the compressor face pressure distortion con- 
tour map. Called by MAINLP , SYMBLE assigns a symbol for each 
discrete point at the engine face, depending on the pressure 
indicated at that point by the interpolation routine, INTERP. 
The spelling of SYMBLE was selected to avoid passible con- 
flicts with library functions in some compilers. 


24, Subroutine TRIDIA 

Subroutine TRIDIA salves the tridiagonal matrix equa- 
tions in the turbulence modelling computations. TRIDIA is 
controlled by subroutine FINITE, which sets up the finite 
difference equations to be solved by TRIDIA. 

25. Subroutine UNSTDY 

This subroutine is the primary subdriver responsible 
for the predictive evaluation of the peak instantaneous dis- 
tortion. Called by the main driver of MAXIDYN, UNSTDY con- 
trols the input and output of dynamic data, manages the com- 
putations leading to the peak distortion prediction, and 
controls the output of results. Some of the more important 
activities and functions of UNSTDY are listed below: 

a. Call NEWPSD to input identification and data control 
parameters for the test run. 

b. Compute compressor face averaged dynamic pressure and 
Mach number, inlet vortex properties, and other para- 
meters leading to the peak distortion prediction. 
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c. Call subroutine LNPQUT to print some of the tables of 
distortion data, and the distortion contour map. 

d. Call subdriver NEWPSD to read in and analyse the dy~ 
dynamic data or select the turbulence modelling rou- 
tines if there is no dynamic data in the input file, 

e. Control the subroutines which iteratively evaluate 
the most probable peak distortion level and print the 
results . 


26. Subroutine TURBUL 

TL1RBUL is the subdriver responsible for the turbulence 
modelling prediction when there is no dynamic data in the 
input file. TURBUL controls subroutines CUBIC, INITL, I NIVEL , 
FINITE, TRIDIA, and PRNT in the synthesis of simulated dy- 
namic data for processing by the subdriver UNSTDY. TURBUL is 
called by subroutine NEWPSD when the user specifies the "no 
dynamic data” option in the input data file. TURBUL assigns 
a finite-element grid to . represent the discrete points on 
the compressor face for the finite difference analytical 
scheme. The boundary conditions for the turbulence model are 
estimated based on the total pressure measurements from the 
steady- state data, and the initial values of the turbulent 
kinetic energy and dissipation rate are found via INITL. The 
inlet face velocities are then found via INIVEL, and the 
turbulent equations are solved by FINITE. These result in 
estimates for the rms total pressure fluctuation levels, 
which are then fed back to NEWPSD for the computation of the 
most probable peak instantaneous distortion. 
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II.B. TYPICAL DATA RUN SEQUENCE DESCRIPTION 


A typical data run of "the MAXIDYN dynamic distortion 
prediction program can be traced as follows. 

The main driver initiates the data input sequence with 
the reading of control parameters and the inlet pressure 
probe configuration. After the radial and angular locations 
of the probe rings and rakes are read in, subroutine NEWPSD 
is called to read in the engine filter frequency, the rms dy- 
namic data cut-off frequency, the data-point dwell time? or 
the length of time in which the dynamic data is measured, 
and the specific distortion factor which the user selects 
as primary for the peak instantaneous distortion analysis. 

The main driver then reads in any data identif ication 
titles which the user may elect to input. The resulting set 
of comments are printed at the top of each page of the out- 
put. If an "end of file" or "endjob" instruction is entered 
at this point, program execution is aborted. After the ti- 
tles are read in, the main driver reads in the base radial 
profile and steady-state pressure array. The steady-state 
pressure recoveries are stored into the instantaneous array 
as a starting point for the peak instantaneous computations. 

Subroutine PRNT is then called to print the table of 
steady-state pressures and the base radial profile. These 
items are included on the first page of the output. Subrou- 
tine PFIX is called next to compute the face-average Mach 
number. Subroutine DISTRT is then called to control the com- 
putation of the steady-state distortion factors. 

Subroutine DISTRT is a subdriver which manages distor- 
tion factor computations. DISTRT computes the simple distor- 
tion parameters before calling subroutines INTERP and DISPAR 
to calculate the eight primary distortion factors. Subroutine 
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PRNT is called if there are an insufficient number of probe 
rings to allow accurate radial interpolation for some distor- 
tion factors, in which case a message to that effect is prin- 
ted in the output. Subroutine INTERP is called to carry out 
circumf erential interpolations at discrete locations on each 
of the probe rings, in preparation for the computation of the 
distortion factors by DISPAR. This subroutine contains a set 
of sample distortion factor which can be modified by the user 
as desired. The results of these computations are eight pri- 
mary steady-state distortion factors used to define the dis- 
tortion level at the compressor face. After computation of 
the steady-state distortion level is completed, DISTRT re- 
turns control to the main driver. 

The main driver then calls subroutine PFIX to set up 
certain parameters required to develop the distortion con- 
tour map for the steady-state case. A "dummy ring TT of probes 
is set up to enable interpolations through the centerbody . 
PFIX then tranforms the input pressure measurements into pa- 
rameters used by the map-generating subdriver routine MAINLP. 
MAINLP controls subroutines INTERP and SYMBLE which generate 
the symbols in the distortion map. INTERP interpolates the 
pressure at the discrete locations on the engine face based 
on the steady— state input data, while SYMBLE assigns a char- 
acter based on the interpolated value at each of the discrete 
paints. After all interpolations and symbol assignments have 
been completed, MAINLP returns program control to the main 
driver. At this point the steady-state distortion has been 
completely defined and the dynamic distortion evaluations are 
commenced with the calling of the main subdriver, UNSTDY , 

UNSTDY controls the subroutines and subdrivers which 
compute the total pressure fluctuation levels which translate 
into so-called "dslta pressures" . These are added to the 
steady-state pressures to produce the dynamic distortion 
level and the most likely peak dynamic distortion. After all 
of the dynamic calculations are completed, UNSTDY returns to 
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the main driver to start another data run. Before this fakes 
place, however, the main subdriver, UNSTDY, manages all of 
the dynamic data input and calculations, or the dynamic data 
simulation if the user selects this option. 

After setting initial values for some of the vortex 
properties, UNSTDY calls on NEWPSD to read in inlet operating 
parameters and some program control parameters . Most of the 
inlet parameters are non-functional, that is they are not 
involved in distortion computations, but rather are of inter- 
est for identification and comparison purposes. One control 
parameter, NTUR, allows the user to select the no-dynamic-data 
option, or to input the required dynamic data conventionally, 
After checking the control parameters for errors, assigning 
default values if necessary, NEWPSD returns to UNSTDY, which 
evaluates the flow velocity and Mach number, and the vortex 
properties. LNPDUT is called to print the simple distortion 
parameters and flow descriptors, and the distortion factors 
for the steady-state case. After printing the steady-state 
vortex properties, UNSTDY again calls LNPOUT to print the 
steady-state distortion contour map. Subroutine NEWPSD is 
then called to begin the evaluation of the most probable 
peak instantaneous distortion. 

After printing out the identification and control pa- 
rameters, NEWPSD branches according to the dynamic data op- 
tion selected by the user. If the no-dynamic-data option has 
been selected, NEWPSD calls the tubulence modelling subdri- 
ver, TURBUL, to compute the rms fluctuation levels which 
would otherwise be input as data. After setting up a fine 
grid at the compressor face for finite element modelling of 
the flow distortion, TURBUL controls subroutines CUBIC, PRNT, 
INIVEL, FINITE and TRIDIA in the development and solution of 
the turbulence equations. Subroutine INITL sets initial val- 
ues for the turbulent kinetic energy, INIVEL computes the 
velocity gradients, and subroutine FINITE uses Finite differ- 
ence formulations to solve the turbulence model. CUBIC per- 
forms cubic spline interpolations and TRIDIA solves tridia- 
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gonal matrices generated by the finite difference equations. 
TURBUL then calculates the r ms total pressure fluctuation 
levels at each of the probe locations, and the results are 
printed. Control is then returned to subdriver NEWPSD , and 
program execution continues as if the dynamic data had been 
input, rather than computed. 

If the user selects the option for the reading in of 
dynamic probe data, then these data are input at this point 
in program execution. In either case, NEWPSD then calls 
subroutine FRF to evaluate the vortex core size as a func- 
tion of the f i ltered- to-unf i ltered rms pressure fluctuation 
level ratio using an iterative scheme. FRF calls subroutine 
SUMMER to evaluate the error function of the inverse solu- 
tion. After the vortex core size is found, the results are 
printed by NEWPSD. When all of the dynamic probe data have 
been read in, subroutine MAXDP is called to compute the 
most probable maximum instantaneous distortion levels. 

MAXDP is the subdriver which controls the computation 
of the most probable maximum peak in the distortion level , 
After determining the mean values for the vortex core size 
and fi ltered- to-unf i ltered rms pressure fluctuation ratio , 
and the effects of engine filters and vortex flux rates via 
subroutines RSIGMA and RATK, the rms and mean instantaneous 
levels are computed by adding a "delta” distortion value to 
the steady-state value. Subroutine EXTRME is then called to 
evaluate the most probable (50% confidence level) extreme 
value of the peak instantaneous distortion. EXTRME utilizes 
subroutine SEARCH which controls the twin subroutines FFXL 
and PFX in the determination of the "delta" distortion used 
to find the maximum instantaneous distortion. EXTRME is also 
used to determine the distortion factor extreme values for 
the 95% and 99.7% confidence levels, using the same plan as 
for the 50% confidence level. After printing the results of 
these computations, program control is returned to NEWPSD , 
and then back to UNSTDY. 
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After re-evaluating -the vortex properties for the max- 
imum instantaneous case, UNSTDY finds the pressure recover- 
ies at each of the probe locations based an the predicted 
maximum instantaneous distortion and the steady-state data. 
Subroutine PRNT is called to print some output, then PFIX 
is called to compute flow velocities and Mach number at each 
of the probe locations. The face-average Mach number for the 
peak instantaneous case is also determined by PFIX, In the 
same manner as with the steady-state case, subroutine DISTRT 
manages the computation of the distortion factors given the 
total pressure recoveries found by UNSTDY for the maximum 
instantaneous distortion case. Subroutine MAINLP controls 
the printing of the pressure distortion contour map far the 
peak instantaneous case as for the steady-state case, and 
LNPOUT is again called to assist with the printing of out- 
put of the peak instantaneous data. After all the output 
for the test run has been printed, UNSTDY returns to the 
main driver. The main driver checks for additional sets of 
data or new test cases. If there data, then program execu- 
tion begins with the reading in of data titles and the 
steady-state pressure data for the new case. If an END OF 
FILE or ENDJOB command is encountered, meaning there is not 
an additional test case. 

In summary, the MAXIDYN dynamic distortion program 
computes the most probable maximum peak instantaneous dis- 
tortion given the steady-state distortion data and limited 
dynamic data. After reading in the steady-state pressure 
recoveries and computing the steady-state distortion fac- 
tors, the steady-state distortion contour map is printed 
along with the distortion data. The average rms pressure 
fluctuation level is then used to determine the most proba- 
ble peak instantaneous distortion. The rms fluctuation data 
may be read in, or simulated by the turbulence modelling 
routines. After the prediction for the most probable peak 
distortion level has been made, a new distortion contour map 





is generated to represent this case. 

The results of the calculations in the MAXIDYN distor- 
tion program are printed on several pages of output. This 
material includes a listing of all input data, the steady- 
state distortion factors and parameters, the properties of 
the convecting vortex used to describe the flow in the in- 
let, the dynamic rms pressure fluctuation data and/or tur- 
bulence modelling data, and the pressure distortion contour 
maps for the steady-state and maximum instantaneous cases. 
Details on the input and output data are provided in their 
respective sections. 
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III. INPUT DATA DESCRIPTION 


The input data are divided into three primary groups. 
The first group defines the inlet pressure probe arrange- 
ment at the measurement plane, and some data control para- 
meters . The second group Includes identification titles and 
the steady-state inlet distortion data. The last group is 
the "dynamic data” - the rms total pressure fluctuation le- 
vels - for each of the probe locations. These data may be 
limited to as few as two probes, or omitted entirely as an 
option to utilize the turbulence modelling dynamic data si- 
mulation capabilities of the program. 

The general arrangement and formatting rules for the 
input data are given in Figures 2 and 3. Further illustra- 
tion on the arrangement can be found in the sample problems 
in Section V. The following is a description of the input 
data items, presented in the order in which they are read 
by the software. Items marked with an asterisk (*) can be 
omitted from the input file without disrupting program exe- 
cution. In this case default values are usually assigned, 
or simulated in the case of the dynamic data when turbulence 
modelling has been selected by the user 

The first group of input data include specifications 
for the Inlet probe configuration, data filter frequecies, 
data point dwell time, and a parameter with which the user 
may select the specific distortion factor used in the gener- 
ation of the distortion contour maps. NR and NP are the to- 
tal number of probe rings and rakes, respectively, and RAD- 
LOC and ANGLDC are the radial and angular locations of the 
probes. KD is the distortion factor key used to select one 
of the eight distortion factors available in the program . 
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The time on point;, or data point dwell time, T, represents 
the duration of time in which the rms pressure fluctuations 
are measured and calculated. FO and FCO are the engine fil- 
ter and rms dynamic data cutoff frequencies respectively. 
Further information can be found for each of these variables 
in the detailed descriptions below: 

NR is an integer corresponding to the number of pressure 
probe rings used in the test run. NR should include 
static pressure rings located at the centerbody hub 
and at the outer radius, even if these are not inclu- 
ded in the i nstrumentati on , so that the distortion 
contour map resembles the engine face geometry. If, 
for example, there are five total pressure probes lo- 
cated along the inlet rakes, NR should be entered as 
seven to account for the static pressure probes. 

NF is an integer correspondi ng to the number of pressure 
probe rakes used in the test run. These rakes are 
generally positioned between the hub and the inside 
surface of the nacelle, and are evenly spaced along 
radii around the centerbody hub. NR represents the 
number of probes along the rakes. 

RADLOC is a one -by- NR array of real numbers correspond- 
ing to the radial locations of the pressure probes 
placed along the probe rakes. RADLDC includes the 
radial location of the centerbody hub, as well as the 
outer radius of the inlet at the nacelle inner sur- 
face. RADLOC may be dimensional, or a dimensionless 
fraction of the outer inlet radius. Units may be ar- 
bitrary in the dimensional case. It is noted, however, 
that the vortex dimensions will be in terms of the 
dimensions of RADLOC. See figure 4. 
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ANGLOC is a one-by-NP array of real numbers which corre- 


spond to the angular locations of the probe rakes. 

The units of ANGLOC are degrees, with the top rake 
being 'zero' and with the angle increasing clockwise, 
as viewed from the front. See Figure 4. 


KD is an integer with which the user selects the distor- 
tion factor of primary interest in the test run. Of 
eight available distortion factors included in the 
program, one is selected for use in generating the 
peak instantaneous distortion contour map which match- 
es and represents the predicted peak distortion level. 
Definitions of the eight distortion factors provided 
in the program are given in Figure 9 . Below is a key 
for use in selecting the desired distortion factor. 
Entering an integer (1 through 8) effects the selec- 
tion of the distortion factor indicated below: 


1: KTHETA (Pratt S Whitney circumferential distortion #1) 

2: KD2 (Pratt S Whitney circumferential distortion #2) 

3: IDG (General Electric circumferential distortion) 

4: IDR (General Electric radial distortion factor) 

5: KRA (Pratt S. Whitney radial distortion factor) 

B: KA2 (Pratt S Whitney combined distortion factor) 

7: DSPR (Delta [loss in] stall pressure ratio) 

8: ID (General Electric combined distortion factor) 

It is noted that are two distinct Pratt S Whitney cir- 
cumferential distortion factors from two distinct def- 
initions (see Figure 9). The combined distortion fac- 
tors are found by combining the circumf erenti al and 
radial distortion factors. In the case of KA2 , the 
circumf erential distortion factor used in the combi- 
nation is KTHETA, The distortion factors represented 
in this program are only examples - the user is free 
to redefine or modify them at will. 
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T 


This is the dynamic data time on point or TT dwell ?T 
time during which the total pressure fluctuation level 
is measured and the root mean square value is deter- 
mined . The units are seconds, and a default value of 
one second is assigned if no value is input. T may 
be omitted if the no-dynamic-data option has been se- 
lected for all test cases. 

FQ 4 This is the engine filter frequency, in Hertz. FO 
is used in the computation of the mean peak instanta- 
neous distortion levels. The purpose of the filter is 
to remove the effect of engine speed on the measured 
pressure fluctuations. A default value of 500 Hz is 
assigned if no value is input. 

FCD ' This is the low pass cutoff filter frequency used 
when measuring the filtered rms total pressure fluc- 
tuation levels. The ratio of f i ltered- to-unf i Itered 
mean square pr essure fluctuations are used to predict 
the most probable maximum instantaneous distortion. 

The units of FCO are Hertz, with a default value of 
1000 Hz, when no value is input directly. 

The second part of the input data includes title 
blocks, the steady-state total pressure recovery array, and 
several inlet flow parameters. TITLE 1 and TITLE2 provide 
space for 160 characters of identifying comments. PS is an 
NR by NP array of steady-state pressure recoveries. The base 
radial profile BRP is the ratio of ring-average pressures 
to the face-average pressure . ALPHA and PSI are the angle 
of attack and sideslip angle respectively, and the f reestream 
Mach number is given by M0. The flow velocity at the engine 
face is U2. BF , CKP , and RKP are weighting factors used in 
the computation of combined radial/ circumferential distortion 
factors. The mass flow ratio, MFR, gives an indication of the 
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the mass flow rate before and after inlet duct bleed-off. 
NTUR is a control parameter which allows the user to select 
the option of inputting the dynamic data, or having these 
data simulated by the turbulence modelling scheme. Finally, 
SPTRC is the total pressure recovery through the inlet shock 
system in a supersonic inlet. 

More detailed descriptions of the data items in the 
second group are given below. Most of these data may be de- 
leted from the input data deck, without causing any real 
difficulties. Many of these are simply included for identi- 
fication purposes, while others are provided with default 
values to avoid data errors. Default values are included in 
the detailed descriptions below: 

TITLE1 and TITLED are alphanumeric hollerith arrays used 
for test run identification. Two lines of up to eighty 
characters each are available for information such as 
engine/inlet type, Mach number, angle of attack, yaw 
angle, altitude, and so forth. TITLE1 and TITLES are 
printed at the top of each page of output for easy 
reference. TITLE1 also is used to check for an END OF 
FILE or ENDJOB command at the end of the data file, in 
which case program execution is stopped. 


EFP 


is the 
defined 
ring to 
array w 


compressor face base radial pr 
as the ratio of the average pr 
the face-average pressure . BRP 
ith a value at each of the radi 
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of the array are NR rows by NP columns, 
are pressures at radial locations RADLOC 
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rows of PS are static pressures assoc i- 
static pressure rings located at the 


28 



centerbady hud and the surf ace of the inlet at the 
engine face. These static pressures can be measured or 
computed values. The pressure array is eventually used 
to generate the distortion contour map, and also is the 
basis for finding the instantaneous pressure array, P. 

ALPH is the a ircraf t/inlet angle of attack relative to 

the freestream, typically in degrees. ALPH is used for 
run identification, and does not enter into any compu- 
tations . 

PSI is the a i rcraf t/inlet sideslip or yaw angle relative 
to the freestream, typically in degrees. Like ALPH, 

PSI is of interest for identification and analysis, 
and does not enter into the computations. 

MO is the freestream Mach number. Of interest for identi- 
fication of test runs, MO does not enter into calcula- 
tions. 

U 2 is the flow velocity in the inlet at the compressor 
face, in feet per second. 

BF is the b-f actor used as a weighting term for the com- 
putation of the combined radial and circumferential 
distortion factor KA2 (Pratt S Whitney). BF is multi- 
plied by the radial contr ibution , and the result added 
to the circumferential distortion to get the combined 
distortion factor. BF has a default value of 1. 

CKP is the circumferential weighting factor used in the 
computation of ID, the General Electric combined ra- 
dia 1/circumf erential distortion factor. CKP is multi- 
plied by the circumf erent ia 1 distortion, then added 
to the radial contr i but i on . CKP default value is 16.4, 


23 



RKP is the radial weighting factor used in the computa- 
tion of ID, the General Electric combined radial/cir- 
cumferential distortion factor. RKP is multiplied by 
the radial distortion factor, and the result added to 
the circumf erential contr ibut ion . The default value 
for RKP has been set at 11.1 in the program . 

MFR is the inlet mass flow reatio, defined in terms of 
the streamtube geometry. Specifically, MFR is the 
ratio of actual inlet mass flow rate, to the maximum 
inlet mass flow rate. The maximum inlet mass flow 
rate is defined as the product of the freestream vel- 
ocity times the inlet hilite area. Low MFR implies 
a large amount of inlet spillage. MFR is generally a 
function of the engine thrust level and flight velocity. 

NTUR is a control parameter which allows the user to se- 
lect the turbulence modelling dynamic data simulation 
capabilities of the program as an alternative to using 
measured dynamic data. Inputting a value of 1,0 for 
NTUR causes the program to branch to the turbulence 
modelling routines within the program. A value of zero 
or defaulting the input of a value causes the program 
to branch to the routines requiring the input of dyna- 
mic data. 

SPTRC Is the total pressure recovery through the inlet 
shock system of a supersonic inlet duct. For subsonic 
and transonic inlets, SPTRC is equal to one. If SPTRC 
is unknown for an arbitrary supersonic inlet, it can 
be estimated by using Figure 5, with a value of 0.90 
being reasonable as a rough preliminary estimate for 
most inlet configurations. A default value of 1.0 has 
been set in the program. It is noted that SPTRC should 
always be less than or equal to one. 
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The third part of the input data consists of the dy- 
namic data. These data consist of rms total pressure fluc- 
tuation levels from fast-response total pressure probes. The 
number of dynamic probe data sets in each run is fully under 
control of the user - within certain limitations. In a nor- 
mal run, the number of dynamic probes is equal to the number 
of steady-state total pressure probes. In a reduced dynamic 
data run, the number of dynamic probes can be anywhere from 
two to as many as would be used in a normal run. Finally, if 
the user selects the no dynamic data option by setting NTUR 
equal to 1.0 [ see previous page), these dynamic data may be 

completely omitted from the input data. 

There are four input variables in the dynamic data. 

NPG is a run identification code, NPR is the probe identifi- 
cation code, RS is the filtered to unfiltered ratio of mean 
square pressure fluctuation levels, and SIG is the root mean 
square level of total pressure fluctuations. Further details 
on these data can be found below. 

NPG is a code number for identifying data runs . This user 
definable integer can be completely arbitrary, though 
entering a value of zero, or defaulting the input sig- 
nals the end of the dynamic data set. Therefore NPG 
can be any integer greater than one. When all of the 
dynamic data has been input, entering a value of zero 
for NPG (or leaving It blank) will signal the program 
to move on the the next phase of computations. 

NPR is the numeral designation for the location of the 
dynamic probe. This identification code can be found 
with the aid of Figure 4, which is given as an example 

far the convenience of the user. Other inlet probe and 
instrumentation configurations may result in a differ- 
ent numeration scheme, so Figure 4 should be used as a 
guide 1 i ne . 
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RS is the ratio of filtered to unfiltered mean square 
pressure fluctuations. RS is found by squaring the 
ratio of the rms pressure fluctuation level filtered 
at the cut-off filter frequency FCO, to the unfiltered 
level. This ratio is used in the prediction of the 
maximum instantaneous distortion level, and a default 
value of 0.50 has been included in the program. In 
addition, a maximum value of 0.70 has been set to a- 
void errors in certain computations. These values are 
easily modified, if necessary, by the user. 

SI G is the unfiltered root mean square value of the total 
pressure fluctuations measured by the dynamic (fast 
response) total pressure probes. Generally, the units 
of SIG are identical to those of the PS array, which 
are nondi mensi ana 1 total pressure recoveries (local 
total pressure divided by the freestream or inlet lip 
total pressure) . 

Sample problems have been included in this manual to 
illustrate the arrangement of the input data, and to further 
clarify the utility of the various capabilities of the pro- 
gram . Figures 2 and 3 show formatting rules and the general 
arrangement scheme of the input data. 
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Figure 2. Batch Input Data Deck Formatting Arrangement 
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Figure 3. InpLjt Data File Formatting Arrangement 

































IV. PRINTED OUTPUT DATA DESCRIPTION 


The printed output of the MAXIDYN dynamic distortion 
program consists of from five to seven pages of data ( sixty 
lines on each page) depending on the options selected by the 
user. The first two pages pertain to the steady-state dis- 
tortion characteristics, with the steady-state pressure ar- 
ray, vortex properties, distortion factors, and related 
parameters, along with the steady-state distortion contour 
map. The next page or two involves the dynamic data, with 
rms pressure fluctuations levels, turbulence data, and 
some additional flow parameters . The following page is a 
listing of a summary of distortion factor extreme values 
as predicted in the Melick prediction technique. The final 
two pages are similar to the first two, but pertain to the 
peak instantaneous distortion level. 

The following is a page-by-page description of the 
output. Since the content of the output depends on the dy- 
namic data input option chosen by the user, some variables 
described may not apply to a specific case. Data affected 
by the dynamic data options are so indicated, and all of the 
affected data are found in the middle pages with the dynamic 
data groups. Some data are provided with default values as 
described in the input section, and the default values are 
repeated here for convenience. 

Page 1 

The first page of output consists of five tables of 
steady-state distortion data. Immediately below the title 
block provided by TITLE1 and TITLES, the steady-state total 
pressure recovery array is printed. The rows in this matrix 


36 



are Identified with the RADLOC radial probe locations, and 
the columns with the ANGLOC angular rake locations. The sta- 
tic pressure rings associated with the innermost and outer- 
most RADLOCs are not included in the pressure array. The 
numbers of the pressure array P are otherwise identical to 
the input array PS. 

Beneath the pressure array is a table of base radial 
profiles. These are defined as the ratio of probe ring aver- 
age pressures to the average pressure over the entire engine 
face. For each radial location RADLOC [including static 
rings) a value of PTR/PTA is given. This BRP array is iden- 
tical to the input array BRP. A default value of 1.0 is 
assigned for each term in BRP when no value is input. 

The next table is a listing of the overall flow des- 
criptors. These simple distortion parameters are used to 
evaluate the distortion factors and vortex properties. The 
terms appearing in this table are defined below: 

PTMIN [also FTMN and TMMI N ) This is the minimum total 
pressure recovery value from the pressure array, P , 
exclusive of the static pressure data. 

PTMAX (also FTMX and TMM AX ) This is the maximum total 
pressure recovery value from the pressure array, P , 
exclusive of the static pressure data. 

PT A VG [also PT A V and TMAVG 3 This is the face-average 
total pressure recovery from the pressure array , P , 
exclusive of the static pressure data. 

PS A VG (also SM A VG ) This is the average value of static 
pressure from the static pressure data in the pres- 
sure array, PS. The two static pressure rings at the 
centerbody hub and outer radius of the inlet supply 
these data . 
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QAVG (also QAV ] This is the Face-average dynamic pressure 


recovery, computed as the difference between the aver- 
age total pressure recovery and the static pressure . 
Mathematically stated, QAVG - PTAVG - PSAVG. 

The three remaining terms in the flow descriptors table 
are algebraic manipulations of PTMAX, PTMIN, and PTAVG, These 
terms are self explanatory - for example, ( PTMX- PTMN ) /PT A V 
is interpreted as the difference between the maximum and 
minimum total pressure recoveries, divided by the average 
va lue . 

Following the overall flow descriptors table is a table 
□f flow distortion factors. The eight distortion factors 
listed in this table are representat i ve of a variety availa- 
ble to the industry, and are intented as examples. The user 
is free to redefine the distortion factors within the pro- 
gram . Next to the distortion factors in the table are some 
weighting factors used in calculating combined distortion 
factors. The eight distortion factors and their weighting 
factors are described below: 

K- THET A (also KTHET A , K6 , and KTT A ) Pratt S Whitney 

circumf erential distortion factor [#1] - see Figure 9. 

KD2 Pratt S Whitney circumf erent i a 1 distortion factor 
[#2] - see Figure 9 . 

( I DC ) -MAX General Electric maximum circumferential dis- 
tortion factor - see Figure 9. 

( I DR 3 -MAX General Electric maximum radial distortion 
factor - see Figure 9 . 

ilf jA Pratt S Whitney radial distortion factor - see Fig- 

ure 9 . 
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K A3 Pratt S Whitney combined radial/circumf erential dis- 
tortion factor - see Figure 9 . 

□SPR Delta [loss in) stall pressure ratio - see Figure 9, 

ID General Electric combined radial/circumferential dis- 
tortion factor - see Figure 9 . 

B-FACTOR [also BF ) Radial weighting factor used in com- 
puting KA2 - see Figure 9 . 

BSF Intermediate weighting factor used in computing ID - 
see Figure 9 . 

KC (also CKP) Circumferential weighting factor used in 
computing ID - see Figure 9 . 

KR (also RKP) Radial weighting factor used in computing 
ID - see Figure 9 . 

Beneath the flow distortion factors table is a 1 ist 
of vortex properties. These properties are described below: 

THMN (Theta Min) Angular location of the probe rake with 
the minimum average total pressure recovery, in deg- 
rees. The 'zero 1 rake is the upper vertical rake. THMN 
is one of the ANGLOC angular locations, and depends on 
the steady-state pressure array. 

DKMN The average pressure recovery along the rake desig- 
nated by THMN. 

THMX (Theta Max) Angular location of the probe rake with 
the maximum average total pressure recovery, in deg- 
rees. See THMN, above. 
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RKMX The average total pressure recovery along the rake 
designated by THMX. 

DTH (Delta Theta] The angular difference between THMX 
and THMN. DTH = THMX - THMN . 

THETA (also THE ] The angular location of the rake mid- 
way between the rakes designated by THMX and THMN , 

THETA = % ( THMX + THMN]. 

A1 (also ART, etc.] The vortex core size which fits 

within the boundaries of the rakes designated by THMN 
and THMX. A1 represents the size of the steady-state 
vortex . 

G 1 (also GAMMA] The orientation angle of the steady-state 
vortex. This is used to satisfy the amplification of 
the steady-state distortion level by the vortex field, 
in determining the peak distortion level. 

Page 2 

The second page of the printed output is the pressure 
distortion contour map for the steady-state case. A repre- 
sentation of the high and low pressure regions at the com- 
pressor face of the enging, this map is useful in identify- 
ing and visualizing the nature of the distortion of the flow 
through the inlet. Symbols are used for identifying the 
pressure at any point in the measurement plane, and a key to 
the mapping symbols is provided. The numbers provided in the 
key are interpreted to mean the percent difference between 
the local pressure and the face-average pressure - for exam- 
ple, an indication of -3.0 is interpreted as three percent 
below the average pressure aver the engine face. 
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Page 3 


The Third page of output; is related to the dynamic da- 
ta, which may be included in the input file, or simulated by 
the turbulence modelling scheme. The content of this page 
(and sometimes the next page) depends on whether the dynamic 
data is input or simulated, as described in the cases below: 

Case 1 : Dynamic Data is Input 

Immediately following the title block is a listing of 
several inlet parameters, along with some parameters used 
with the dynamic data. These parameters are described below: 

T^ The dynamic data time-on-point, or dwell time during 
which the dynamic data are measured for each of the 
dynamic probes. The units are seconds. 

FO The engine filter frequency, in Hertz. FO is often 
associated with the engine rpm speed. 

R_T The outer radius of the inlet at the compressor face, 
or the location of the outermost static pressure 
probes - the maximum value of the RADLDC array. 

HI The inner radius of the compressor face, the radius 
of the centerbody hub, or the minimum value of the 
RADLDC array. 

ALFH (also ALPHA) The inlet angle of attack, relative to 
the freestream, in degrees. 

PSI The sideslip or yaw angle of the inlet in degrees. 

SFTRC The total pressure recovery through the inlet shock 
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SPTRC (cont’d) system in a supersonic inlet configura- 
tion. For supersonic inlets SPTRC is less than 1.0, 
while subsonic and transonic inlets will have SPTRC 
equal to 1.0. 

MO The freestream Mach number. 

ETA The face-average total pressure recovery from the 
steady-state pressure array, PS. 

MFR The mass flow ratio of the inlet system. This gives 
an indication of how much of the inlet air remains 
after bleed-air has been removed, 

U2 The inlet flow velocity at the engine face in feet 
per second. 

QPT2 The dynamic pressure divided by the total pressure . 
QPT2 = QAVG/PTAVG, where QAVG is the face-average dy- 
namic pressure and PTAVG is the face-average total 
pressure , (see "overall flow descriptors" table des- 
cription in Page 1 descriptions.) 

RS AT FC = The cutoff frequency of the rms dynamic data, 
(see FCO in input data descr i pt i ons ) 

The next table of data includes the dynamic data as 
input by the user. For each dynamic probe location selected, 
values for the rms pressure fluctuations and the resulting 
vortex core size are given. The specific terms in this table 
are described below: 

PROBE The numerical designation for the location of a 
dynamic probe. See Figure 4 
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RS The ratio of the filtered to unfiltered mean square 
total pressure fluctuation level. 

SG /PTE (also SIG ) The unfiltered rms total pressure fluc- 
tuation level. 

A/RT (also ART ) The mean vortex core size, based on the 

magnitude of SIG. The vortex core size is nondimension- 
alized to the inlet radius, RT . 

Immediately below the dynamic data table, the average 
value for the rms unfiltered total pressure fluctuation le- 
vel, SIG, is printed along with the average vortex size. 

These terms are actually used in the Melick peak distortion 
prediction technique (Reference 5). 

In some cases, specifically when the number of dynamic 
probes in the dynamic data is relatively few, the distortion 
factor extreme value table is printed on page three immedi- 
ately below the dynamic distortion table. The reader should 
refer to Page 4 output descriptions for identification of 
the terms in this table. 

The following is a descripti.on of the terms on page 3 
of the output when the no dynamic data option is selected 
in the input file, that is, when the dynamic data is simula- 
ted by the turbulence modelling techniques: 

Case 3: Dynamic Data is Simulated 

The data appearing on the third page of output includes 
all of the data appearing in Case 1, excluding the dynamic data 
listing. The reader should refer to the descriptions in Case 
1, except for the dynamic data - PRDBE , RS , SIG, and ART. 

These terms are replaced by three tables of turbulence cal- 
culations, and the simulated values of SIG for each of the 
available dynamic probe locations. 
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The Following is a descriptive listing of the data 
on the third page of output when the dynamic data are simu- 
lated by the turbulence modelling technique. The parameters 
listed in Case 1 are included in these data, and the reader 
should refer to the description listing there for details. 

The first table following the inlet and control para- 
meter listing gives dimensionless velocities of the flow at 
the compressor face for each of the pressure probe locations. 
These velocities are calculated based on the steady-state 
pressure data from the input file. The rows of the velocity 
array are associated with the ANGLOC angular locations of the 
rakes, while the columns are associated with the RADLQC ra- 
dial probe locations along the rakes. The first and last 
columns reflect the static pressure probes located at the 
centerbody hub and outer inlet radius. 

Following the table of dimensionless velocities for 
each of the probe locations is a listing showing the itera- 
tion of the turbulent kinefic energy, and the kinetic energy 
dissipation rate. The relative error in these terms is mini- 
mised during the iterations. The first column shows the 
error in the turbulent kinetic energy; the second, the error 
in the turbulent kinetic energy dissipation rate, and the 
third gives the sum of these two . These errors should de- 
crease rapidly within the thirty iterations allowed. Once 
the errors have been minimized, the turbulent kinetic energy 
and dissipation rates are used to generate the synthisi zed 
dynamic data . 

The next table is a listing of the results of the tur- 
bulence calculations, including the synthesized dynamic data 
and the turbulence modelling parameters. This table is simi- 
lar to the dynamic data table in the Case 1 descriptions, 
but includes some additional terms. This table is large, and 
may actually be slipped to the fourth page. The terms appear- 
ing in this table are defined on the following page. 
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PROBE The numerical designation of a dynamic probe, used 
to define the location of the probe. See Figure 4. 

UU The sum of the squares of the radial and circumferen- 

2 

tial velocity gradiants in (ft/sec) 


E The turbulent kinetic energy dissipation rate in units 
of f t^/ sec^ . 

2 

K The turbulent kinetic energy in (ft/sec) . 

SG/PT2 (or SIG) The synthesized unfiltered rms pressure 
fluctuation level . 

Printed below the synthesized dynamic data table are 
values for the face-average SIG, and the mean vortex size, 
ART. These values are used to predict the most probable 
peak instantaneous dynamic distortion level. 


Page 4 

The fourth page of output is a table of distortion 
factors and parameters leading to the most probable peak 
instantaneous distortion for each of the eight sample dis- 
tortion factors. The terms and distortion factors appearing 
in this table are described below: 

KTT A (or KTHET A ) The Pratt S Whitney circumferential 
distortion factor, definition #1 (see Fig. 9 ). 

The Pratt S Whitney circumf erentia 1 distortion fac- 
tor, definition #2 (see Fig. 9 ). 

1 PC The General Electric circumf erent 1 a 1 distortion fsc- 
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-tor (see Fig. 73- 

IDR The General Electric radial distortion factor (see 
Fig. 7 3 . 

KRA The Pratt S Whitney radial distortion factor (see 
Fig. 73 . 

KA5 The Pratt S Whitney combined radial/circumferential 
distortion factor (see Fig. 73- 

DSPR The loss in stall pressure ratio (see Fig. 93. 

ID The General Electric combined radia 1/circumf erent ial 
distortion factor (see Fig. 73 

STEADY STATE This column indicates the steady-state 
values for the distortion factors, as computed from 
the input distortion data. 

MEAN VALUE The mean instantaneous distortion level, com- 
puted by adding the mean instantaneous rms fluctuation 
level to the steady-state distortion. 

SIGMA INF The unfiltered rms distortion fluctuation le- 
vel . 

SIGMA FO The rms distortion fluctuation level, filtered 
at the engine filter frequency, FO . 

MOST PROB The most probable peak instantaneous distor- 
tion level at a 50% confidence level. Statistically, 
this is the most likely value for the peak distortion 
level in the statistical prediction analysis. Moving 
away from this value decreases the probability. 
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confidence level percent: of area under curve to left 


95% PROB The peak instantaneous distortion level at a 
95% confidence level. This is interpreted as meaning 
there is a 95% chance the actual peak instantaneous 
distortion level will be less than the indicated le- 
vel. The likelyhood that the actual peak will reach 
this level is small. 

99,7% PROB The peak instantaneous distortion level at a 
99.7% confidence level - there is a 99.7% chance that 
the actual peak will be less than this level. It is 
very unlikely that the actual peak distortion level 
will ever be this high. 


The most probable peak instantaneous distortion level 
far the distortion factor selected by the user is printed 
immediately below the distortion factor extreme value table. 
This distortion factor is used to develop the peak instan- 
taneous pressure array and contour map. 


Page 5 

The fifth page of output is much the same as the first 
page, except the data applies to the peak instantaneous dis- 
tortion rather than the steady-state. The terms in the tables 
are defined in the descriptions of Page 1 , though any refer- 
ences to the steady-state case are understood to be replaced 
by the peak instantaneous case. 

A major difference between the fifth page and the first 
page is that the vortex properties table has been replaced 
by a vortex location table, with some new terms. These are 
defined on the following page. Many of the terms are similar 
to some of the steady-state vortex properties, though they 
apply to the peak instantaneous vortex. 
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VBAR The average vortex strength in terms of the vortex 
tangential velocity vector nondimensionalized by divi- 
ding by the Flow velocity at the engine Face. This 
property is used to deFine the source oF the pressure 
Fluctuations . 

A/RT (also ART] The vortex core size in terms oF the 

vortex radius divided by the inlet radius. See AY/RT . 

GAMMA A vortex orientation angle in degrees, due to the 
rotation oF the vortex core about the x axis. See 
F i gure 1 0 . 

BETA A vortex orientation angle in degrees, due to the 
rotation oF the vortex core about the z axis. See 
F i gure 1 0 - 

AY/RT (also AYRT) The vortex core size in terms of the 
vortex radius (the radius at the maximum tangential 
velocity oF the vortex system) divided by the inlet 
radius. AY/RT is also the same as A/RT. 

VL The nondimensional vortex length limit. VL has been 
set at 999.999 (inFinity For all intents and purposes) 
in the current program, though this is easily altered. 
VL should represent the true vortex length limit divi- 
ded by the inlet radius. 

ft/RT (also RRT) The radial location oF the vortex cen- 
tral core. R/RT is a dimensionless value with a maxi- 
mum value oF unity. 

THET A The angular location oF the vortex center in de- 
grees. Zero degrees is the top vertical position, with 
positive THETA being clockwise about the engine Face. 
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VBMAX The maximum vortex strength in terms of the tan- 
gential velocity of the vortex divided by the flow 
velocity at the engine face. VBMAX also appears as 
VBM in the FORTRAN coding. 

VBO The vortex strength as approximated from the total 
pressure rms fluctuation dynamic data 

AO/RT The vortex core size computed from the average of 
the dynamic data power spectral density (PSD) funct- 
t ions . 


Page 6 

The sixth and last page of the printed output consists 
of the pressure distortion contour map for the peak instan- 
taneous case. The terms and parameters appearing with this 
map are identical to those in the steady-state map. These 
parameters are described in the Page 2 description in this 
section . 


SUMMARY OF DEFAULT VALUES 


The following is a summary list of default values for the 
input/output variables which have such values: 


T 

= 1.000 

BRP = 

1 .000 

F0 

= 500 Hz 

MFR = 

1 . 000 

FCO 

= 1000 Hz 

SPTRC = 

1 . 000 

BF 

it 

o 

a 

o 

VL = 

999 . 999 

RS 

= 0.500 

max RS = 

0 . 700 

CKP 

[or KC) = 15.4 

RKP (or 

KR) = 11. 
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SAMPLE PROBLEM 


V . 
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SAMPLE PROBLEM 


V . 


A. Introduction 

Four sample data sets are provided to illustrate the 
input/output capabilities of the MAXIDYN distortion program* 
These problems are taken from provisional data, and represent 
a variety of inlet operating conditions. The first and second 
cases are supersonic inlets with a full set of 40 dynamic 
probes, and a partial set of 14, respectively. The third case 
is a subsonic inlet with the minimum number of high-response 
dynamic probes - 2. The final case is a transonic inlet with 
no dynamic data Input. This case makes use of the turbulence 
modelling capabilities of the program, which simulates the 
dynamic data. 

Some of the primary data parameters are shown in the 
table of part B, below. Figure 7 gives a complete listing 
of the input data files for the sample problems. Figure 8 in 
part C shows the output from the four sample problems. Defi- 
nitions of each of the terms in the input and output listings 
may be found in the input and output data descriptions of 
section II, parts B and C, respecti ve ly . 

B. Sample Problem Input 

The four test cases provided here have similar probe 
ring/rake configurations. Figure 4 illustrates the arrange- 
ment of the pressure probes at the engine face. Some of the' 
main parameters in the input data are tabulated on the fol- 
lowing page, with a complete input data listing in Fig. 7 S 
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PARAMETER 

CASE 1 

CASE 3 

CASE 3 

CASE 4 

NR 

7 

7 

7 

7 

NP 

a 

S 

8 

8 

RI 

1 .000 

1 .430 

0.383 

1 .645 

RT 

5.169 | 

4.346 

1 .000 

5.000 

KO 

3 ( IOC ) 

3 ( IDC ) 

3 (IDC) 

6 (KA3) 

T (sec) 

30.0 

3.0 

3.0 

3.0 

FO (Hz) 

500 . 0 

500.0 

1000.0 

1000 . 0 

FCO (Hz) 

1000.0 

400.0 

500.0 

500.0 

ALPH (deg) 

4.0 

5.0 

- 

- 

PSI (deg) 

0.0 

0.0 

- 


MO 

1 . 36 

3.5 

- 

- 

U3 (fps) 

514.3 

378.5 

433.4 

575.8 

BF 

0.733 

0.733 

0.733 

0.784 

CKP/KC 

1.0 

1 .0 

- 

! 

RKP/KR 

1 .0 

1.0 

- 

- 

MFR 

- 

1 .0 

- 

- 

NTUR 

- 

0 

- 

1 

SPTRC 

- 

1 .0 

- 

- 

probes 

40 

14 

3 

0 


[A dash indicates a defaulted entry, or zero,, The 

program assigns default values in these cases.] 


C. Sample Problem Output 

The line-printer generated output for the four test 
cases is presented in Figure ©. The number of pages of output 
varies with the dynamic data content, but never exceeds seven 
pages, unless there are run-time errors (exponential over- 
flows, negative square root radicals, etc.) or compile-time 
errors. Run-time errors can occur with bad data. The terms 
given in the output are defined in Section II, part C. 
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Sample Problem 
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Supplemental equations 
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q « average dynamic pressure at engine 
face 

M = number of rakes 
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Figure 9 . Distortion Factor Definitions S 
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vortex orientation angle between y axis 
and the x T -y T plane 
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vortex orientation angle between x f and x 
axes, with the x axis in the x T -y T plane 


Figure 10. Definition of Vortex Angles 


82 


APPENDIX A. 


PROGRAM SOURCE CODE LISTING (FORTRAN) 


83 



im$L PAGE'® 
poor QyMnn/. 


X co o j 

X X » 3 

l/.l m. ft ft 

>— » ft K-< Z! LL 

«•' UO Cl *-• X 

Q x x ^ g x 

LL - X* X G - 

■’• CT‘« - X ft- X 

CO •“■. 4M ft il 

o O ». OJ 21 

~ uj cu > x x 

w x» - x o h- 

X’ X X 1“ — 

> X *-* X - IJD- 
X X * w - ft* > 

X — i/) x X 

U X .iXX » X 

— CI3 — — ft- h* 

ft~ X - - 

i i O I jO ft- 

1-4 X •*-* G *• o 

*■-** X *•-'* X G OJ 

CO - X X ft- 

X X - »x OJ 

O UO — . ft- X 

X 'X« i— • CO x 

ft4 — ' 1""4 I* ft— 

X ! » 7 '« OJ '*-** ►—* i — t 

S*. CO M I* j [«M 

IL v h Q X — 

W Ui J H a 

- X* w ' C5 CO O 

OJ XXX - OJ 

XX - X UD ^ 

K - X - > — 

LU X X X X 

X 1 ** - X X 

— CL OJ ly) 

G G ft- - *-h i 

— ' CO - X IX 1 ft- 

X */.’ 21 X’ *—■* — 

Q ^ > V O 

o - x x — 

OJ ,-•. X X O W 

O H- - ^ G 

O ^ G X ‘ i— » 

G X' G X G 

X X X X X* ft- 

'X G - X CO X 


OJ C IL 

n y. a 

iii □ 

• 3 J Z 

2 J C 2 D 

Q Cl w q h 

M ft M ft 

ft ~ c ft «r 

C I.»’I X C U 

CL *£• IJ Q 

j 2 II □ J 

I Ii h Q J 


Hi 

j 

G 

X 

A 



K 

H4 UJ 

H- 

u. 


ft 

J 

IA. 

ft 

ft 

u 

UJ 

X 

L- 4 

X 

X 



H- CO 

• - X 


ft 

ft 

a: 

ft 

£ 

V. 

r. 

- 

•s 

® 

X 

IX _J UJ UJ 1=1 

Q 

□ 

ft 

44 

3 

□ 

ft 

M4 


5* 


-1— I 

1 - 

h- X 

i 2 ; 

u 



ft 

i£i 

ft 


X 

o 

X 



X 

CO X 

X X G 


CL 

“i 

X 


Vi 

.*• 


-•PH 

X 

OJ 

X 

X 

B 

hhD 

UJ 

ft 

z 

ft 

S 

44 

n 






H- 

X 11 

• • H4 

ft 

ft 



ft 

>1 

G 

G 

OJ 

X* 

I£. 


X X 

ll II CO 

□ 

tr 

ft 

Ml 

i.'i 

z 

M 

G 



X 

X 

CO 

00 

X 

X' X 

HM 1— 

~ a g 

K H- X 

CL 

CL 

ft’ 

X 

ft 

44 

4-» 

J- 

z 

G 

G 


iV'i 

H- 

1-4 

X G 

X X G 




u 

IJ 

ft 

HI 

G 

x* 

.1 

X 

X 

t— 

G X» 

H* H~ hi 

ft 

Ml 

Mi 

□ 

c 

Hi 

-J 

X* 






X X 

CO CO ft* 


44 

44 

ft 

ft 

jr 

UJ 



•si» 


•si* 

•S'J 



Cl 



ft 

jr, 

CL 










(L 

ft 

>1 

2 

*- 

Cl 









H« 

Z 

z 

ft 

C 


,41, 



V 


ft 

X X 



UJ 


X 

G G 



ft 


!J 

9, ». 






H m-4 



z 


(4. 

It tl 



□ 


Mi 

— 1 4*4 



44 


ft 




ft 



k A 



ft 


ft 




□ 


z 

G* *“* 


1— t 

ft 


r 

’w* '•«*■ 


• 

Ml 


□ 

X s X* 


iT* 

44 


ij 

a a 



ft 



XXX 





ft 

G L4 l£‘ 

! 


ft 



- X G 

X 

X 

“i 


C 

XXX OJ 


G 

LL 

.-V 


G || 

II 



1-4 

in 

y-* 

OJ 

x 

44 


ft 

.4^, A-.. "“J 

“i 

a 




' 1-4 -r-4 


X 

a 



® O O 

,*s 

4-4 

ft 

i/‘« 


i-4 1-4 m 4 t— 4 

a 

X 


X 

ft 

- - - ft* 

i— 

X 

ft 

-*ts. 


UO If) UT5 X 

X 

11 

Ml 

X 

c 

ft- 

ft- 

X 

LL 

G il 

ft 

CO X 

CO 

a 

2 

jpM* 

*1 

fi m C h '•-■ G 


ft* 

ft 

uL 

ft 

X X X 11 11 


r i 

2 

-i X 


iii X 111 H IL OJ Ii 

X 


X CO 

ft 

X X X "0 »— • “> 

•— 

X 

J 

X 4— 





ft 


LL 




ii 


2 




1 j 


4H 





1-4 



u u u u u u u u u i ; u i ‘ 1 1 i . 

x ;i* £' H © 'J? n* ‘X 1 0 ° *- 4 oj co -r uo «xt r- o:« -5 —» cu co rf uo *x* r-- co >t. ® — oj co 

!r; ;z; |z‘ ’x jz; ® ® , ccrrccr.r - oj cu oj oj oj oj oj oj oj oj oo co co co 


Q4 



Oii41 READ <5* 1£4»END=999.> TITLE1 

0042 20 PER [i (5» lc'4> TITLE2 

0043c 

0044c INPUT f.RSE RRD I ML PROFILE 

0o45c 

0046 PERU <5*101> < ERP <J> » J=1 , HR> 


ORIGINAL iFASE fS 
OF POOR QUALITY 


Ml 







£ 


% 




LU 


it 


<f-4 


H 




w " 


Ui 




cj 




:> 


a 


s 


£ 


x 

« 

£ 


CL 


0 

•> 

£ 


a: 


x 

c 

LL 


c 


LL. 

EL 





t 

X £ 

z 


a 



x <r 

□ 


H 


N 

Vk 

M 


z 



■y—i jjJ 

H 


1-4 


CJ 

It £ 

£ 




□ 

“ i 3 

□ 


U 1 


-j 

r. 10 

H 


u 


U4 

Ml 

Mi 


£ 


'X 

“ Ui 

w 


J 


Cti 

£ 

A 


CL 


V. 

X LL 


o 



\ X 

X 

Ui 

• 

A 


ZI 

04 Ui 

H* 

T— ! 

Z 



H* 

C 

It 

£ 


’*-« v-4 

- £ 

J 




II 

H 

3 

“3 

U1 


CJ ►“* 

X Mi 

ij 

* w ‘ 

£ 


a 

X 

J 

CL 

3 


— j »• 

«. V 

•r 

X 

Mi 


w 

** A 

ij 

Cl 

Ml 


X H-. 

ll C 



Ui 


CL ~ 

“? Ui 

i 

'2: 

£ 


1 p 

- K 

ij 

O 

LL 



Mi 

w ■ 

■ 



i*/) 

"3 

X X 

o 

UJ 


X 

H 

Z X 

a 

H 


CJ 

H- Z 

X 


'I 


6 

»--« CL m 

Ml t— 

^ LU 

K 

uc 

X 

~ X £ 

CU Uj • 

a 

M'l 


Q ^ X 

Cl 

Z H- 

w. 


e-. 

x o *. 

•w 

^ X 

f — 1 

> 

’I-* 

x ^ ^ 


H* H* X h- 

II •- 

A 

II 

*-• ii 

x r-- h 

X 3 H-4 |V'l 

” X 

c 


II lf> w 

it 

X G X *—* 

CL 

LU 



o H 

X £ X l— * 

U 0 CL 

K 

® 


*- X « z 

A 

cu 

Mi 

CO 

w Q 00 

~ h- — £ 

X I» XX 




h“ X 

~ *-* II CL 

X Mi X X 

Cj Ll 

K 

G LL LU C 

X X 

'X »x X 

Cl ►— 

“! 

Q 

i — i CL G 

X 3 X J 

w J CJ «x 


LL 



J 

J 





<r 

c 


M 



ij 

IJ 

UO 






CU 




CO 


i i 

lj u 



ij ij 1 1 

IJ IJ IJ 

n. *:c 

CU* ^ * 

C»J 

co rr uo 

N. CO CP* C X* 

CU CO T 10 * J f--- 

tt 1 ; 

U“ in U0 

uo to it* 

UO ITt IT' lO ' 0 *•£» 

X* X ; '•£' X «.£i ‘.J 


w w 


; z* 


•*, O o 


' 1 1 

® 

1 ■ Z 




• 



Ml 


Mi 



Ui 


£ 

Ui 


”“[ 

h 

ui 

z 


J 

3 

H 

ll 


£ 

Cl 

Ui 

J 



h 

Z 

Cl 




»r 



3 

□ 

£ 

H 


□ 


<1 

Z 


Hf 

ui 

Cl 

Ui 



I 


z 


£ 

H 

IX 1 

Ui 


D 



£ 


K 

Ml 

HI 

“i 


Ml 

J 

Cl 

Ml 


M 

G 

ll 

'I 


A 

£ 

*r 

1U 



h 

z 

T~ 


it‘i 






5 

□ 

H 


s 

ij 

H 

“i 


Ui 



CL 


3 

0 

M 

z 


£ 

Mi 


M 


H 

J 

Ml 



3 

£ 

h 

*T 


£ 



M 


H* 

;> 

Ui 



Mi 

A 


w 


z 

h 

Ui 

Mi 


H» 

Mi 

£ 

LU 



Z 

in 

6 




*r 

£ 


3 

» 

iu 

Cl 


H4 

hi 

z 




LL 


Z 


£ 

£ 

H 

a 


z 

3 

X ( 

H3 




LL 

H 


LU 

Ui 

z 

£ 


H 

£ 

M 

J 


3 

3 


□ 


LL 

Mi 

z 

ll 



Ml 

£ 

£ 


G 

Ui 

□ 

- iu 


ij 

£ 

lL 

LU H 



LL 

Ifj 

iVj 3 


D 


z 

X 


h* 

M'i 

>r 

X 




£ 

X Ml 


Ui 

a 

H 

- J 

X 

3 

LU 


a 

X 

M 

3 

□ 

X £ 

T” 

h 

£ 

H 

*— * H 



A 


X z 

X Oi 

5 



X 0 

z: n 

£ 

£ 

M 

1 l 

j— • 

A 

H 

IL 

X 

X X 

Z ! 

lu 

LL 

X LL 

x X 

1,0 

2 


X J 

X CO 



J 

CJ 3 

CJ H- 4 

J 


J 




A 

«r 

£ 


£ 

3 

ij 



ij 

>1 

u u 

u u u 

1 l 

SJ 

ij 

iT. <Z 

CU tj* 

UO j N 

CO 

iT. 

X 

N- N N. r j r 

N* N. X 

(X 

IX 

O Z : 

■Z Z Z Z 1 

Z i Z ! Z' 

»z 

Z 1 


85 



CALL UN ST DY 

6n to 10 



i ii 

* A 

✓v — >X 

U~ iJL ® 

»-» 00 cu 

CU 



iT< O «>—4 rj* 

O; O O CU 

u " u 

— « oj co to '£• r-- co *x. 

CO CO CO CO CO CO CO CO CO 


Lt 

>Ti 


86 



009 £ SUBROUTINE PFXL <flKS» DKS» RHT> 

0 093 c 

0094c THIS SUBROUTINE SOMES ''S I <3Mfi K FOR KBflR/SIGMfl K.L.E 


o 

* 

CU 


>— 1 o 



r- 


A 

o o 



o 


X 

cu cu 



r- 


X 




» 

cu 

* 

a a 





A 



i/9 

♦ 

a 



i>*;» 

v’ 


H- 

p 

□ a ■.£ 

X 

X 

i>9 


* 

IX* X* CO 

X 


v 

a 

p 

* 

\ 

rC 

X 


1 


LO 

i—* 



*w* 

o b * 

'■ft 

CO 

co 



■ • 4- 

CO 


in 

o 

*;.* 

1/9 

|V| 

• 

• 

• 

Ui 

• • J*. X 

. 

o 

o 

CD 

1 

UJ LU H- X 

o 

+ 

+ 

# 

O 

j jr ♦ 

-1- 

iT* 

in 

UJ 

a 

■ -'XT r- 


Cu 


<*-* 

I— '>•) '4& T 

jT, 

*-£‘ a 

« 


T m v' X U”> 

o 

\D 

IP X 

•jf 

4* 

>r x c *j- 

o 

• 

* 1 

A 

X 

-*■ '*-* -j • 

« 

O 

w X 

V 

11 

X CU 

o 

11 

11 If 



X Ul i! 11 

11 

g 

X X 

IL 


*—* •— • X 'X P 

x 

X P 

U* 

p 


*.£« r*- QO Ch o 

iT- iT. }T, iT» o 


X 

P Q 
uj x 

Osi IU 


CU 

tj » ! IJ 

cu co rr to «.£» r**. oo cr*» o ^ 


87 



5 c THIS SUBROUTINE CONTROLS THE OUTPUT 

bC SOME OUTPUT IS CONTROLLED BY ' CliHI'PI' ' AND FFI X f . 

7 C HP? RABLDCj fiNGLDCJ HMD P HAUE FIRST ARRAY DIMENSION EQUAL TO 11 . TEN OF 


• 

in 

ft ft 

H X 

Z 

X OO 

p-t 

X wft 

V» 

I/O »v 

Vi 

P -4 . 9. 

Ui 

-• in x 

IJ 

a so x 

0 

X -• cj 

£ 

w. iT‘i 

Cl 



® ,® 

J 

M X OJ 

•I - 

X 

<L. 

i£i CL 

£ 

>'X — 

UJ 

xx y 

1 

£:•>■ UO- 

Z 

O X X 

M 

9. Q ft' 


,-* H- 

£ 

o x O 

D 

^ 'X ^ 

IL 

X ~ 


CO m. CL 

M 

X .*-*. X 

UJ 

□ uo ~ 

M 

X iD ~. 

UJ 

M O 

UJ 

X ;'T‘ CU 

z 

m tyy •. 


X . ‘ ^ ftM 

Vl 

Cm LU ft-' 


I 

h 


Ui 

Vi 


£ 

a 

Ll 

LU 

J 

w 

•I 

J 

M 

c 

i 

Ui 

c 

»r 

UJ 

Vi 

Ui 

X 

w 


IJ Ij 

CO cr- 


cu 


^ £ Cl 

—4 M •/ 

£ x 

CL £ -- 
CL H- X 
» ». »X 
- CU £ 
> CL £ 
'Iflh 
£ I — *• 
X - LD 

•w —4 X 

~ X X 
CL A £ 
** H- H- 
IX w* 
LD K 
£ »■ O 
•X □ OJ 
CL h— 

•>. iX CU 
.vs H* UJ 

i>*;» i 

H- 


O H" h~ 

ax - 

-J h— .*-*.■ 

* , o ^ a oo o 

oj x x x ~ OJ 

X LL «• LC X 
b- t »*. CL •» O *— 
aJ « I ix >I LU 
Cm .O ». - X X 


3 !»•/) 
X 0>0 

V. 

C 

'X 1 j 

CU;~ 

a v 
-J x 


uo 

'•£> 

LM 

X 

X 


OJ 


£ 

C 


l£t 

C 

Cl 


o 

cu 


•X* 


U“ uo 

ft— ft— 

ii !! 


X X 
UJ fj 

X X 
X X 


X 

CU 

I/O 

1- 


uo 

X 


X 



\ 

X 


•ft— 



ft— 






<® 

OJ 


h- 

ft. 

J — 4 


’•D 

X 


X 




l£i 


II 



II 






ft— 

ft— 


X 

t£i 

1— 


*w 

H- 


K 



0 

T" 


►— 



p— 








£ 

0 

h— 

ft 


X 

£ 


p— 




X 





ft*. 








£ 


*>0 

..•S, 


x 

X 


K 



£ 

X 


.ft 



ft. 







■w 



*■- 

o 


<r 

X 





LL 



p— 



— 








X 

X 


ft— 


X 

ft. 





I 

4 











X 

X 

X 

I j 

O 




a 


X 



— 



H- 



H- ID 





H* 

H“ 

1— 


ft-< 

£ 



K 

; 

X 



Cl 

p— 





i 






p— 

p— 

£ 

X 


p— 


iT, 

X 

CU H* 



s --* 

! 


a 



a 

X 





X 

X 

X 

_l 

X 

£ 


ft— 

K 


p— 




X* 


X 



X 

T* 





■*^ 5 : 


X 

X 

l£i 

>— 


**— * 

i>0 


H* 



£ 

•*— 





•«* 

K 







ftv 

X 

I/O 

X 


K- 

9 * 

OJ 




□ 

■- 













X 


ft. 

H* 


G 

H— 

«p. 




u. 

L— 



... 

... 

... 

_ 

_ 






X 

X 


0V 


a 

H 


® 

uo 



X N 

i® 


f*- 


CO 

t® 


uo 

X 


.***. 


ft* 

C* 

X 


X 

X 00 

iT. 

o 


Z ft -4 

a 

1® 

-1— 

o 

»® 

■—4 

CU 

OJ 

co 

1^, 


Z* 

O 

’:*** 

,.-V 

— «* 

H- 



1- 


ft— 



Q 

X 

ft— 

•»— 

ft— 

rf 

•r— 

ft— 

ft— 


ft— 

-X 

•XI 

Aj 

UJ 

o 

^ =r 

X 


£ 

I/O 



r. 


M —4 

X 






m. 

ft*. 

ft— 

ii 

X 

a 


X 

OJ 

X. 

X 


G 


OJ 

•jC‘ 

•ill 


H 11 

ii 

•vD 


*£» 


•■jj 

'•jj 

«.£( 

u 

-*® s 

11 

X 

q; 

*1 

■>4 

o 

H- 


p— 

X 

1 




IJ — » 


■•ftp** 

'w* 

• w * 

*■— 

w 


•w* 

p— 




X 

OJ 

x» 

t£i 

ft 


iV'i 

X 





Ui 

p— 









Z' 

H 

» 

& 

1 

X 

i/O 

m 


£ 

x* 

X 

X 

X 


Vi OJ 


X 

X 

X 

X 

X 

X 

X 

z 

1® 


— * 


ft- 

X 

ft 

X 


X 

p— 

*w 

H 

1— 


Oj K 

1— 

I— 

h~ 

H- 

1— 

H- 

H- 

ft— 

ft— 

X 



ft— 

i/.i 

H* 

L -4 



!£» 


p— 

p— 


h 

t 

H— 

p— 

p— 

p— 

•— 

p— 

p— 



ijji 




X 

X 

J- 


—ft 

G X X 

X 


z a a x 

X 

X 

X 

X 

X 

X C G 

X 

X 

Ll. 






A 

X 

p— 

X 

■»• 


LL Pi 

X 

z 

y' 

3 

3 

3 

3. 

3 : 

♦— i 

iZi 

X 

(ft— 

►— 

vV 

ci* 


;-v 








h 

























0 











Z‘ 











c 




cu 










Cip^ 


*z 









OJ 




cu 










ft— 


ft— 





u 






u 

Ij u 














OJ 

m 

rt 

UO 

*•£ 

r*. 

co 

t T. 

o 

ft— 

cu 

oo tT li” 

r L‘ 


iT 

>T. 

■Z 

ft - 4 

OJ 

ro 

rr in 


x 

CO 

Oj 

OJ 

CU 

CU 

CU 

cu 

cu 

CU 

CO 

CO 

co 

; * t i • « i * i 

CO 

iT 

oo 

.V i 


TT 

*T 

rr 


ft— 



T 

— 4 

ft— 

ft— 

—4 

ft— 

»— 

ft— 

■— i 

— *< 

ft— 

ft— 

ft— ft— -r— 

ft— 

ft— 

ft— 

ft— 

ft— 

«— 

ft— 

ft— 

ft.. 

ft— 

ft— 

ft— 

ft— 

;® 

o 

>Z 

•z* 

o 

X- 

I - — 1 

o 

c* 

o 

I*™*- 

■Z l O Cm 

o 

w 

«®. 

o 

z 

Z 

z* 

t® 

. ® 

»® 

w 

z 

Z- 


cu 


X 


X 


fA 

X 


X* 

X 


1.U 

X 


«T‘ 

cu 

X X 
G X 


UO uo 


88 



OUTPUT SECTION FOP USER-DEFINED I'lSTORTION PFIRftMETERi 


ORIGINAL PAGE IS 
OF POOR QUALITY 


\ \ \ 

rr rr 




nr - 





X 

a a a 

>r 



X 








X 

0 - P- P- 

IL 











X 1 

XXX 












i >0 

a. 4 . bn 

— r 











00 

11 11 11 

□ 











X 

X > > 

P 4 











X 

X X X 

H* 



j- 








X 

1 — 1 — ! — 

K 












XXX 

O 

^ uo 

to yo -5 

.-/I 








X 

\ \ \ 

K 

a s 

a * 49 ^ 









l£t 

( p., /»., 

Ui 

4-4 --4 

«-4 4 - 4 %^ 









X 

r 1 r 

M 

4-4 4-4 

4 P 4 4 - 4 -^. 




i— *■ 





X 

xxx . 

A 

X X X x €2 









X 

hhh P-- 


fk as 





4 * 






X X X X 

M 











X 

ill- 

Ui 











X 

X X > 11 

z 

V 



p- 







'■£' 

X X X X 

HI 

X 

X 




+ 





4-4 

\ H“ L— H* >* 

Ll 




o 







V. 

X X X X X 

Ui 

1 

- •» 

\ 

LU 


II 





X 

t/\ *-■’ •*-■■ 0 

M 

\ 


3 

X 







0 J 

X X X X X 


CO X 

X- • 


»y;i 



4-4 


Jj-. 


CO 

a p.. p-. p-. uo 

Z 

X CU O 0 J il® 


Cl 



a 




X 

“ 4-4 ^4 —4 C*. 

Ui 

□ LU •— 1 

X 




il 

l D 




X 

X. 9, w. 9. *:•:* 

Ml 

1 — **-*• 

X LU ^ 


BP, 


iVf \ 

X 


+ 


1-4 

^ X X X IT* 

3 

>•_' XXX 

\ 

X 


J 1 

w. 


X 


X 

cl 0*1 ro »V, ^4 


|T* 4-4 4-4 

--4 4 - 4 -^^ 


lu CL 

a 

X 

\ 

IT.' 


□ 

Q a 

IL 

X <44 4-4 

4-4 



-L 

LQ 

< 1-4 

</> 



X 

iV'» 00 00 I*V| 1*0 

□ 

•s 9. 

- r-X: 


r» 

X 

X 


X 

X 


X 

LU • • • • 


XXX 

X X:t^ 



«* 

> 1 

*— i 

X X 



G P*- P-- pw 

1 - 

L 1 4-4 4-4 


*44 

>.£‘ 

X 

»> ’* X 

0 J 

X 




LL X X X 

3 

H 4 4—4 4*4 



V* 

X 



'-D 

f -4 


IJD 

3 «• K k a> 

il 

H* \ \ 

\ \ < 

X 

—-4 

CJ 

X X 

in 

X 

X 



- □ u it ii 11 

H 

X UTj to 

to to-'^ 


II 


X P- 


G 

1-4 


X 

X X X X I-D 

3 

a - *• 

. .m 


H 4 

X 

»-4 uo 

•X 


uo 


X 

X X ^ X > > 

□ 

1- ^ ^ 




l /. ; 

X 

X 

i“r 



\ 

0 X X X X 


> 4-* 4-4 

<o -4 «a* 4 -*rt? 

yc 

o*. 

Lu 

X e- 

H- X 


"x* 

,/4 


CU II — i — H* </} 

Ui 

H 4 ll O ll IJL 

\ 


4. 

>X X 

4ZH .A 

0 

X 

G 

x« 

*• X X X X X 

I 

9, 9. 


\ 

H 4 


X 

X •£> 

’•£' 


CO 

•"* 

X X X X X X 

H 



X 


CO 

■*-4 

0 :=* 

CO 

+ 

<—4 

X 

\D X l -£« «.£‘ '•£> '■£' 





£.1 

w 

Q ^'..-4 

X X 




X 

\ UJ 9> »w «•. r. 

£ 

O X 


If} 

ill 

Uj i 

X \ 

X X 




1 

"T > X X X X 

O 

X X 

ctx 



X \ 

> O-a 

LU X. 

x x a- 

X X- 

u* cu x a 

CU \ 

X ~ ♦ 

rf x X 
*"« N- CU 

■w co - 

hhh \ 
X X X , /■ 


X I + i — 1 + *— » 

*- • I N I LL 

x ll '£> »£' OJ 

i-h '•-** ll ^ X »* X 

X X ,< “ t **■ ‘-L 1 

tJ“ <—4 «s iT. fy 

— - CP* X CO 

v. x 1***1 *x* \ x 

X X OJ \ OJ 

oj cu oj \ co 


X Q tO U” UO IO 
*=' X OJ OJ oj OJ 
CU TP 


’L> J<„ 
w- CO 

*-* o*i 

X \ 

CO -*-» \ 


X X X 


r-- 

X l— H* 
- X X 


Ll X X Q fr 

X h~ 1 fr-* - 

CO LU o O 

OJ X X s X 

CO 

X I w ct | g » ®>. 

^ X *-• X lQ ^ 00 X 

- I X X X X IL IO 

\ it—* 4-4 4-1 <*-4 ««. OJ 

\ 4-4 4-4 4-4 4-4 X \ 

* W «*- "* *> w ’ 

X X X x . . X 


LL lL lL. 

a a c- 

X Ll ll i>-;» 


LL LL LL 

OCX 
Ll Ll Ll 


X Cl X 

a o 

LL LL 


tT £■ r X *T* x ^ O ^ cu oo 

; X O 4-4 4-4 |“U |“U Cu CU 


rf uo \0 r r CO *T’ X ^4 Oj CO TT U”' '-£• P-- CO CO G -r-« OJ CO rf U0 
U U” U -■ U - 1 U” ll” '•£ '.£t ».£i <£i '.£ *.£* «.£» ».£* '•£* ’.£ P*- p*- P*- P- P- p*- 


1 L 

t— 

4-4 4-4 

44 iVb^BT ^ -S. 

X 



X 

4-4 -4 

X H“ 

X 



Lm 


X X X 



□ 

X 


0£ U’ 1 X 

H* 

u«<* 

M 

c 


O c 

LU 

X 

H 

X 


%, U. t~ 

X 

X 

IJ 


v|i 

--i> .^1 -r i* ; i< 



Ui 






I.P, 

£ 


a 6 

•£' 


ij ij 

co 


^ *T 

1 ; 


P-- CO 

CP* 

•^T' 4-4 

oj <r u‘< '•£. r-- 

CO *J< 

c 

r- r- 

P-- 

CO CO 

i*Ci 0 * ffi sT: rn 

CO CO 

iT. 


O O G O O G <p G» 


89 



I 

I 

I 

I 


I 

\ 

i 




I 

I 

I 

i 




i 

i 




— “l 


— 



s 


rt 



X 





LL lj 


□ 



»■ 2, 


1 — 



LL, >~* 





Ll H* 


a 



LL “» 


l£i 



LL □ 





X 


£. 



LL E 


■ 



H* Ui 





X H 


• 



Cl j 





Ht 


LU 



UJ U. 





X 


a 



>— • Lu 

x 




H- X 

• 

CL *•£» 



X H 

C- 

x r- 



□ 

II 

u. - 



X Ml 

iLL 

‘w , 



Lu »_i 

X 

*-4 



X 1 

JL 

X ti 



«>*> Ml 

Cl 

' Cm i 



i-* 





I 





h 



IJ 

IJ 

lj ij »J 



0 j 

jv* 

Lf ; \C f*-- 

iV 

iT: X*, 


iT. 

iT, iT, iT. iT. 

tfr. 

tT. X . 



•*— <4 i—4 1—4 

X 

«**■» r»j , 



O o Q 'X* 

X 1 



□ 


a 

in 


® Ui 

LL 

u:» x 

'■£' LL) 

- I 

K O 

in . 

a £ 

— ' LLj h- 
o lm LL 

• *x »>♦> o 

— U. « X 

H ♦ X II 

x A v jc. 

tL 1 U- 


LL 

P Q 

LU X 
LL UJ 


^r 

ij 

*T if* •.£» r*-- 


ij ij 
CO »T, 

•\j <:u 


90 



ORIGINAL PAGE IS 
OF POOR QUALITY 


H 

I J 


LU 

□ 


u. 

Cl* 

LU 


CL 

LU 


LU 


LL 


E 

□ 

E 

E 


Ui 






Ui 






X 




t — r 


H 




• 


Ml 




CO 

m 


LU 




II 


J 

O 

CO 


UJ 


□ 

CO 





Ml 


■X' 


■ 


LU 

K 

CO 

•>. o 


05 

co 

LU 

♦ 


X 

•X’ CO 




M 



+ 


O 

h 

X 

£ 

II 


X h- 

\ i£> 

1 


E 

mi 

Mi 

w 

I 

H 


LU 'X 
X H 
~ X 


2 

OJ * C 
t'lj <V, 

ll a ii ii 

^ S Z Ui 


LU 


i! 


a 

in 


hi 

OJ 




'•u 

a 


c 

in 


Jl^ 

X 


•X* 

a 


a 

x 


ij 

X UJ 

CL X 

LU X X 1 
X* 


II 


w CU ui 2 
LL H- J| X I 


® -J O 
• « • 

5C **— * 
i II 

: ll 

.* LL Cl 
:> m lu 


r- 

□ w 


{_.*• 1 
4> 


♦ CO 

*”* CO 
Cl T. 

LU i T> 

♦ 

X* lx 
CO OJ CL 

:o i uj 

LU I 


_~S 8«i?5w2£ 

°j 2 ft ♦-iiuzs . x . c 

c ^ □ ii. n r u. s ii 

'X' m 'X s i— ^ ir *- . r- i= 


ii 


Cl Ll H»»s ^«»4 • ■,• • 

lL 


ii 

LU H- if - |{ 


ll 

Ui X 

ir 

X LL LL 

II U_ a LL. UL Ui. 

h~ 

Lu 


*.•** ■ 

>x *-* u LU t— 1 LJ 

X 


P O Q 

U 1.1 u It "' J T ~ ~ •-• 

-X 2 2 i !2 i U o: ' *' ?; r <u « rr U-. ■.* r-. CO ,* .= - ,-u u - , a , ,, , T . . ' L ' ' ' " u '- T 

^888888SiSSg|S|S|||g§§||||||||||||||||||| 


T\ 

*r ip 

Oj OJ 


91 



ORIGINAL PAGE IS 
OF POOR QUALITY 


LL »x* oj 

X P » X' 

I/O ft. ..ft,, ft. ft. 

*“ l ». ft -4 X X 

ft. \j*j 1*1 H X 

a ‘X* *x x ll. 

LL *■• UL X »♦ 

•>. tj*» ^ it t- x 

/■v co p *. x 

O 1 O ». OJ x 


Mi 





■h UJ CU 

> X zz 




K 





fj «>. 

X P H 




*m 





IJJ CX 

X H” * 




h* 




tt 

> .X 

X * X 




0 




in 

X LL 

ft. ft"4 C-“ 




X 




Z 

y *. t'/'i 

H-* X X 








M 

O Cl Cl 

X O X 


s'**. 


P 




Ml 

- a - 

•*. 1— K 

a 

’it 


HI 




Ml 

h- 

UJ *♦ r. 

Ui 

Ui 


X 




Ui 

O Q : — 

I£| h- ..-v 

H 

H 


l£l 




lj 

^ X — 

X ^ 

*“i 

n 






Q 

LL Y 

x a cu 

D 

a 


Ui 




X 

CO »* tX 

X h- 


*■-*■ 


X 




Cl 

X X 

r. X OJ 

It 

a: 


H 





a m - 

1- X 

o 

i”i 






J 

x x 

— </} X 

• 

• 

W 

X 

H 

Ui 

Ui 

'I 

P * O 

*■4 «■> H* 

T-* 


a 

(fl 

Mi 

X 

X 

Z 

q iT’t n j 

■*-* ^4 ►— 1 

l 


w 

“! 

it 

K 

1- 

X 

co ft 

<J h- h- 


0 

H 

□ 

M 



Ui 

X '**' *+ 

a x ^ 

B 

\y 

<1 

It 

IL 

Z 


H 

P X — 

X 1- 

Iv 

a 

ij 

X 


□ 

• 

T 

ft. CJ ’*-* 

g !>■:• o 

a 

it 


H 

111 


Ui 

M 

OJ X X 

X - OJ 

it 

LL 

Ui 


X 

X 

0) 


XX* 

X IJD w ’ 

it 

LL 

*>» 

4i» 

M'l 

t- 

111 

3 

z 

H- * X 

ft -»“4 

it 

C 

Hi 

UJ 

B 

H 


□ 

X zz 

XXX 

- »r 


H 


ill Hi It Z 

Z 2 □ 

C *1 « Z P 

^ LJ 


P ‘ 1 ■* •> X -J 

^ il C<J MO y 
Q *'-** X H- ** *-* 


1 1 


D 

It 

P 


Cfc 

a 

cj 

‘X 


Ui 

X 

H 

a 

H 


Ui 

K 

M’i 

Ml 

'I 

in 

U 

X 

c 

1 

X 


ij lj 

ft-4 1*1 j 

UO Lf» 


u 

CO rf 

in in 

nj i*i j 


U U 
U“ ’X 1 
uo u“ 

CO CU 


U i j u u n 

r-- co on 'X' ^ 

uo ir> u* \D 'X* 

CU Cu Cu Co cu 


u u 

CO Cm "T LO ‘X 1 r- CO CT % ‘X 

•X' x x x x x x x r- 

Cu CU CU CU CU CU CU CU CU 


I j ij ij u ij i ij 

ft-* f!_i CM rf U~ ‘X 1 N CO l 7‘* O CiJ 0* 

r*- f s - r>- r-~ r-- r*- r*- r*-- r- co o:> co co 

CU Cu CU CU CU Cu CU Cu CU Mj cu Cu CU 


'■ii 

cu 


x 

*x» 

UJ 


i 

CO 


u~ 


H 

x 

y 

UJ 



M 

u. 


^4 

Ui 

M 

X >'/J 

y 

* 3 * 

H -4 


X 

*1 



3 



X 

’•L 






j 

UI 

'Z ; 




It 


-l 

Ui 

r. ^ 

y 

ft. 

i/.i 

."•S, 

ft -4 

H 

:> 


t?i 



a 



z 




c 

X 

• 


I 

h 


D 

it 

UJ 

■’2 ^ 

r. 



X 


U 1 


l 




H- 

co 


cu 




z 

K 

X 


IJ 


M 

H 

»r 


B*. 

X X 


i 

»>•;» 

I 

i 

1 

Z 



•X 

I 






D 


; n 


HI 

Ui 

3 


j 


Oj 

X 

» 3 ‘ 

O 

*•-*’ 

X 1 

h 

i 

*v‘ 

□ 



X 

ft -4 


X 




H* 

z 

• 


I 

z 

•I 

J 

H 4 

Ml 


K 

ft. 


zz 

X 


V 





LL 



y 

X 



H 

ft# 


X 

yr 

in 


•t 

c 

M 

C i i—i 

y 

X 


HI 

X 

p 



iT. 



❖ 




\ 



Z 


■ 

X 


z 

o>. 

3 

*■•. 


□ 

»x 

y 

X 

y 

ft. 

«T 


Z 

co 






X 

c- 



UJ 

p 

y 

y 

Ui 

3 

IJ 

Z 

•r 

z 

X X 

X 

X 

13 

l— 

\ 

c 

z 

Ui 





■Z' 


l£; 

**-* 



2 

ua 

X 

H 4 

2 

M'i 

□ 

Ui 


□ 

13 y 


X 

»y;t 

X 



UI 

H 

Ui 



X 

Cu 

»'*■ 


i/i 



Z 

lit 

X 

X 

M 

* 

J 


Ui 

HI 

T* 

p 



h* 

X 

z 

H 

3 

z 

iT*> 

ft -4 

ft -4 

1 

4 - 


P 



(— 

w 

3 < 



h 

UJ 

l£. 

-T 

z 

H 

X p 

u 

X 

.-V 

ft. 

ft -4 


*”i 

□ 

HI 

CO 

uo 

1 

a 

> 

O 

X 


co 

IJ 


X 


3 



£ 

»r 

HI 

«p. e. 

H. 


o 

X 

\ t 

z 

□ 


J 

ft. 

ft. 

o 

.i 

+ 

0 



X 


Z 

» 

X 

a 


*r 

M 


Ml 





K 

X 

□ 


Z 


ft -4 

ft- ' 

OJ 


x 

o 

X 


r-ft 

z 

z 



n 

H 


M} 

Mi 

□ 

X o 

LLI 

o 

w' 

X 

X 


z 


z 

u 

II 


y 

♦ 

c 

X 


uc> 

0 

H 

• 

f— I 

Li 



4 n» 

■^r 

It 


■•/ 

OJ 

X 

X 

X 

z 

o 

• 

z 

c a h* 

X 

x 


y r*- 

0 

z 

M'i 

X 

r- 

3 

£ 

ij 

Ui 

0 





1—4 

K 

1 

□ 

■ 

i— » 

hi 

H— 

y 

X X 


X 

X cu 


z 

> 

3 

y 

Ms 

Jm 

□ 

z 

Hi 

t£» 

X X 

CU 

o 

in 


p 

HI 

ft-4 

1 




a 

X 

y 

a 

II 


II 

IJ 

i.i t 

• 

a 


M 

j 

HI 

K 

z 

c ^ 

y 

X 

i/\i 

CO 


h 




O 

X 

X 

X 

X 


X c 

X 


X 

X 

1£* 


L» 

A 

Hl 

HI 



X 


X 

X 

c 

• 

4 

. • 
l\ 

44 

— * 

LL 


t y 

w 

in 

y 

\n 

z 

X 


X 

Z 

UJ 

*r 


M'i 


y y 

T— 1 

i/i 

h— 

p 

1 — 

J 

x 


z 



■*-* 


;/[ 


Ui 


X 

w 

K 


i 

HI 

X 

z 

*> 

□ 

HI 

c o 

y i 

X 

X 

y 

X 

""t 

o 

Q 

□ 

□ □ 

It 

11 

II 

X 

X 

c 

X 

IJ 

Z 

X 

i/.i 

V 

H 


!r 

it 

«r 

ijjl 





p 

U 

z 

Z 

h 

HH 

P 


X 

X 

H 

X 


a 


X 

ft*"4 


HI 


r. 

z 


z 

cc 

•:U 

•:£ 




J 

it 

LL 

1 1 










*L 

h 



z 

z 

It 

z y 

Ui 







•I 

LL 

Z 

»E 










X 

1 



w 

HI 

z 

»E 


Z 







u 

>t 

it 

Z 










u 

Z 



IJ lj lj 

X' fv CO 


rf X 

co CO CO co 
CU CU CU Oj 


92 



‘X X 
X X* X' 
X* UJ X 

x 5 5 


® o UJ I- ♦ 

_• * - X X X 

w O O «. • H 

• » ■ X X UJ 

LU I— H— h* H* X 

X -1 X UJ LU h- 

* • • X X t! 

X X X h- I - X 

w *-*' '•-* '■-* K 

LU 

X UL X X X X 


□ 

I 

IL 

«r 


Ui 


• 7 *' 

O 


a 

l + 




a 

+ 

co 

IT, iT. 

CO CO 




X 


cu cu 




,A, 

r~ 

■*— J 





uo 

N. f“- 




X 

a 

- * 00 






rr 




*w‘ 

11 

X 

II II CO 
X -X OJ 




o 

a 

h- 

h* H— CU 




X 

XXX « 




© 

XXX x» 




X 

f— 

» — 1 — 11 




X 


X 

© 



• 


b- 

■ 

00 


H* 

© 

o : = ( X 




IJ0 

• 

- - X 

II 

CO 


• 

© 

® o K 

X 

00 


X 


a tt 

K 

cu 



i— 

H- x 

X 

• 


X 

X 

X X* CO 

X 

X' 


a 


8 8 

1 — 


Ui 

a 

o— 

CO 


o 

I/I 


« 

a c lYl 


■ 

J 


Cx 

L «4 l. u 

*-* 

* 

Ui 


X 

X X X* 

o 

x» 

co 

It 

(J 

a 

• 

• • • 

ft 


Q 

X 


X 

X 


LL 


h 

n 

□ 

It 

w 

i?i 

J 

0 

w 


W 

w 

C 


H 

>1 

J 

□ 

iL 

It 

Ui 

H 


X M 

»"* J 

J 
£ 

I 1 


X 

X I 
X 
X 

X 

II 

X X 
Ct J □ 

x x n 
) — x 
x > a 


a 

x x x H- 

X X o 

x x x a 

«X O X X 


X 

X 

X 

II 

/“w 

b 

a 

M 

X 


© - LU 


w i— . LU 


X 

o 


r— M 
X X 
X X 


Ij ij IJ 

cu co *t uo X' r-- co o** 

»T’ iT. (T*. tT. IT. |T. fT. jT. iT, 

>‘U CU CU OJ CU CU cu cu cu 


IJ IJ IJ 

Cu CO UO 


iT. 

X 1 r-- co 


IJ 

iT. c* 


■H Uj 

|V ( ( *r, 


IJ IJ 
CO rf 


93 



CO 

4^ 


II 

0 

fl 

II 

0 

0 

0 

0 

0 

u 

j O' 

5 u 
0: 
0 

0:; 

O' 

O' 

O' 

: o: 

« 0: 
0 
0 
0 
0 
0 
0 
0 
II 

i ; : 

ii: 

i 0 
0 
0: 
fi: 
0 
0 
0 
0 
0 
0 
II 


315c 
3 1 6 
3 17c 

31 Sc 
J 1 9 c 
<z Or 

32 t c 

-•c idC 
■« c ! C 

IJ 3 £4 


>3 0 c 
S 3 1 


334c 

335c 

336 

"•3'dC 
3 3 9 c 
34 0 
341 

34 2 

343 

344 
>45 
346 
347c 
;l4'dC 
3 4 9 c 

35 0 
35 1 c 

5 c C 
353c 
354 


SUBROUTINE HI SI FT 

SUEPOIJT INF. WHICH CALCULATES SIMPLE D I STDRT I DN PARAMETERS. 

MF< PAR OC) ANGLDC > AMI' F HAL'E THE FIRST 

ARRAY IiIMEMSIDM EFL'AL TO 11. TEN OF TH E 

11 POSITIONS ARE ANA I LABLE FOR USE. THE 

REMAINIMS POSITION IS NEEDED FOR INTERNAL PROCESSING. . 

COMMOH AB»AN6LI1C<S0) » AUD» BETAS* BF*CDPDP3<10> » CMAVG C10> jFDj ishft » 
ii: CUM I n 1 0> »BKMX<10> »DKSS(10> »FACE<39*65> > FACTOR* FARCE <39* G5> » BSF > 
:l J l , J£ , KEY > KB * LCHTFL * MH * HP > HP » P <1 1 * £ 0> * PR < 1 0> > PS < 1 1 * £ Co » CKP > 
ii-. PSPEC '£0> , RVALUE » QAV » QPT2 » R * PAUL DC < 1 1 > » RAHbE >RI> RMAV > BRP <11 > » 
ii PT * SGDK < 1 CO t SGP < 1 0> j S 1 6 r SMAVG » ST AT I * ST ATO * T * TDP 1 * TDF'£ > TMM I H » 
ii, TDP3 * THETA » THTM I H < 1 0> > TITLE1 (SCO > TITLES <20> * TMAV6 * TMMAX > RKP > U£ 


D I MEHS I OH SEG C 1 0> > RADP <8 j l CO , THTTOT < 1 0> * 
LOGICAL STAT I * STATO? LCHTFL » LADD 


:mhax < 1 0> 


CAL CUL FIT I ON OF SIMPLE CIRCUMFERENTIAL DISTORTION PARAMETERS. 

DO 10 I = J 1 » J2 

CRLCULRTIDW DF CMMflX?CMMINJ HMD CMRI/G FOR ERCH RING 

CMtTAX C I > = 0 . 0 
CMM IN < I > =50000. 0 
CMAVG < I> = 0. 0 
DO 11 J= 1 » HP 

CM MAX < I > = AMAX 1 <P < I » J> > COM AX < I > > 

CMH 1 H < I > = AM I H 1 CP < I j .0 > CMM 1 H < I ) > 

1 1 CMAVG ( 1 > =CMAVG <I> + CP C I > J> .' FLOAT CHP> > 

CALCULATIDN OF CDFOF'S FOR EACH RING. 

1 0 CDP0P3 < 1 > = CCMAVG I > -CMM I H C I > > /CMAVG C I > 

CALCULATION OF TMMAX AND TMM I N . 


tmma::=o. o 


ORIGINAL PAGE IS 

OF POOR QUALITY 


■ 

l£i 

£ 


X 

















I { 














cu 
















£ 














a 
















Ui 



























































£ 












o 


a 
















c 

















X 












u. 












C 1 


X ! 



of 

N-4 
























x» 





y 

y 









X 

in in 














1/.! 





21 

y 









X 

V 

V 


M 












4- 





1— 

l— 










X X 


z 










X 


».u 





e. 

tr. 










y y 


p 












X 

y 

* 













y y y 


X 




CO 






l 


II 

cu 

’’■S' 



*— 4 

t— t 









\ 

\ \ 


H 










i.U 



J’ 



•w‘ 

•w' 






__ 










a 




X 


x 


X 


X 




X 

»n 





y 


li*1 

21 

X X 


Ui 




H- 




X 


II 



a 

X 



•X 

H-4 






ST, 


• s 

H-t 

M *- 4 


\r 








II 


i.U 



y 

1 




y 

<1 







> Il 

X 

X X 


•r 




□ 




a 


X 


X 

X 

8 



T 

y 

X 







>-« p 

X 

y y 


j 




X' 




X 




X 

y 

o 


rj 

t“i 

h 



CU 


X 


P h 

y 

y y 


D 


o 










* 

X i 

X 

• 

cu 


*w 




y 

CL 

+ 



1 

1 1 


ij 

0J 

fH 




■ 

X 



X 

i.U 

H- 


* 

o 

— *t 

<^■4 

f-4 

Ul 



V. 



o 

X It 

x 

*X* IJj 


J 

y 

r- 

o 

<r-4 


X 

X 


X 

X 

*x 

X 1 



o 

«■- 


X 

C 




9. 

y 

■ 

x« c 

X X > 


•I 


■*—4 

■ 

• 


CO 




• 

' w ' 

£ 


H—i 

o 

s 

X 

t— » 



■**"> 

y 

•*— < 

CO 


*x» 

X 

X X 


ij 

•*— » 

11 

o 

X* 


X 

r-4 

X’ 

f ‘t 

H * 

X ! 


* 

X 

p 

y 


X 

X 


• 

if 

II 

y 

+ 

y Z 

y 




y 

X 

II 

ill 


X 


X 

c 

Xi 

a 


y 


UO 

II 

x 

X 

c 


X 1 

H-4 

y 

X 


X □ 

y 

y K 


□ 

ll 

X 


• 


X 

J 

8 

X 

• 

X 

+ 

X ! 


il 


11 


H* 

B 

II 



il 

►— * 

II M 


*-• 


H 

♦— i 


X 

X 


» 

4 

X 

\n 

— ' 0J 

1 1*1 


■ 






H 

— i 


CO 

rn 



X H 

11 

It II 




CO 

X 

X 


II 

CU J_ 

X 

X 

■4“ X 

X 



X 

*—* 

N 

X 

»— 

*1 

II 

•y 

0‘ 

CO 

X 

II 

> ir 


CU CO 



iT. 

-*-4 

*•••* 



1— 



y x 

y w 

X 

*'**' 



X 


y 

y 

— i 


Jy> 



iVt 


X J 

XXX 


□ 



i n 



1-4 

II 


u u 

l! 

II 




y 

g 

y 

y 

_■ 

t— H 

y O C K 

M 

y n 

v— * 

p p 


M 

c a x x 

It 

x a a x 

—-4 t— 1 

CU X 

cu 

X 

X 

X 

y 

Pi 

y 

y 

ij 

P 

X 

p 

P 

X P 

l — U 

y 

1— y 


h 


t— 

j/;. 

1—4 

X 

X 

u- i X 

H— 1 

x X 

X 1— 

X 


i«— 






_j 











»J 



















X 







»r 




Ui 



















IJ 






CO 

u 




Ml 



0*i 















r**- 

IJ U ij 






CO 

ij ij 

i i 


IJ 

U 

ij 














IT* 

\D 

N 

CO 


CU 

CO 

y uo 


CO x ® 

«*— * »CL1 

co y 

in 

••£« 

0- CO 

iT. 

O 


cu 

co 

’T Ul 

*•£' 

r- co 

iT, P> 

•— * 

cu 

C 1 ** 1 

y 

UO 

U“* 

UO 

uo- 

th '■£ *.x 

*0 

•X' 

'•X 

•■L* 

'•L> 


X* X N 

r- r- 

r- r- 

r- 

N 

r-- N 

N 

co 

iTi 

CO 

i*0 

CO CO 

cc* 

co CO 

CO «T‘ 

iT. 

\T. 

iT, 

'7' 




{Vi 

r”» r ” j"t“ 








|*M 


t Vi 

co 

/V, -**| 


Tm 

|V; 

co 

*y 

ft* ,■'* 


CO m"i 

CO C"* 

co 

rn 

co 

i* “i 

C 1 1 

o 

P 




O 



o 


. O P 


o o 

o 


O 1 


p 

o 


o 


® 

® ® 

® '= ! 

P' 

p 

p' 

p 


95 



!£• IS 


\£\ 

OC 

o 


OJ 


rf ** X 
^ h h 

'X G LL 11* 


IlQX>I 
»— CD •—» X 


CU — 
OJ CU 

u~ x* r- co 

tT. iT. iT. j-. 


ir, 

2 

a 


□ 

w 


D 

»-* 

H 

s 

H4 

tt 

K 

z 

D 

!J 

l/t 


•w ••-•' w* w i r, w* | r, 

\n |£, |jT| in Ijr, 

> > > > a: > x 

x x gc x x x x 

O ijj fji ijj q w 

V O I ■ & * & 

hOKh ,-v I- 
JUJ J J ID 

■ o « • 

/■». /“s. *-* CU 


... w- ... a. w * x 

x x x x i x i 

• ■ • « 


H 

IJ 

U1 

Hi 



I 

_4 

^4 


2 



IJ 

X 

X 

X 

IJ 


CU 

‘I 

1 

1 

1 

•r 

rn 

OJ 

Ui 

CU 

CU 

OJ 

Ui 

OJ 



X 

X 

X OJ 



a 

2 

11 

It 

11 X XX 

u. 

X 

H- 

G 

X 5 

o 

o »>. *1. 

a 

K 

' 


X 

X 

X >-* j— * 



x* 

2 

I—* 


1^4 -■• 

VJ 

a 

!£; 

t-* 

X 

X 

R 

p 

p 

2 

in 


2 



1 1 

a 



W 

.*4.. 



w 

,a. 


2 


/**. 

ym, , M ►— 4 

H 


i— 4 

H 

1—4 

t— t 

t— 1 t— « HH 

Z* 

1—4 


2 

K 

Z 

D 

U 

U 1 
2 
H 

2 

Ui 

l 

H 

UI 


in 

X 

o 

• 

H- 

in 

« 

a-* 

X 


rr 

OJ 


Ql ©* *. 


CJ *• */') CU *■ o 


rt C* 


X *-* 11 ~ 


y 

w 

2 

X 

y 

c 

Ui 


Ui 

€ 

-X 

Ui 


€ ill 
2 

X 


a 

j 

Ui 


l2i 

a 


©* 

A 

i=» O OJ 

A 

H 

Q 

iX 


a 

cu 





X 



X 

in 

i- 

*X 

X 

XXX 

X 

X 

H 

-j— 


X 





.***•. 

X 


z 

X 


dim 

X 

X 

X X - 

X 

«K 


o X 






*f 


»■/» 


□ 

B 


Ui 

• 

« 

s a >— * 

a 


P 

CU • 



a 



l/.l 


+ 


1*4 



2 


-y 

w 



A 




H- 



X 

in 



Mi 

t— ♦ 

>— 4 

tn 

>— « 

»-4 

H H4 U, 

1— 1 


«r 

O t-4 


t— t 




it 

X 



5/1 

*•-" 

*w 

ill 




•w 



H- 



a 



l/.l 

i/» 




IX* 

iX» 

Mi 

x 

X 1 

in i£i \ 

X* 


a 

X* 


X ( ^ 

OJ X* 



X 

4 s 

X s 


ct 

J* 





> > 



w 

X > 


> OJ 

4-4 





X 


iu 

X 

X 

I 

X 

X 

X X -» 

X 

4—4 


X 1 X 

o 

X 








Ui 

X 

X 

ij 



XXX 


X 


T* 

X 

X c 

a o 



•2' 


11 


M 

i’«j 

o 

i 

Q 

CJ 

O *X 1 

X 

1 

□ 

Q 

t— 4 

O h- 

H* • 



<9 

H* 


1 


« 

0 

Ui 

B 

B 

a • r»j 

• 

cu 

2 

•2‘ a 

X 

e 

O 



o 

C 


</'• 

h 

Q* 

X 


H 

K 

>3 K X 

K 

X 


• H- 

+ 

x a 

a • 



« 

h"- 

■w 

X 

M* 

X 

l£t 

mi 

X 

X 

x x 

X 

'w' 

Ml 

o X 


IX* X 

?X« »3 

*3' 


X 

H* 

X 

11 

Ui 

« 

B 

UJ 

• 

» 

a a V 

a 


Ui 

a a 

X 

8 

X 

« 

o 

X 

X X 

{■/) 

ifi 


» 



.“S 

- H 


11 

2 

*3* 





4—4 

e 


(V* 

-'.xl 

i 

4— * 

a x 

M 

CU 

CU 

OJ Cu O 

OJ 

o 

M 

X OJ ix 

a U 

j-* 

H 

0S 


11 



»r 

X 

x X 1 

2 

X 

X 

XXX 

X 

X 

2 

a 

X 

x • 

• -i 


4—1 

X 


c- 

•Z 

j 


h 

>1 

X 

H 

X 

C 

X 


u 

IJ 

i 

w 


tT. 


Ui 

• 



>-4 


X 


X 


Ui 


X *-* 11 4 - 

X 


CU X x« 


2 l 



X 

t n 

h- i 

I© 

tt 


w 







s.. 



H 


►— * .-.. ■*-*' 

■ 


* 

X 


X 


K 


X 

h- X 

X 

D 

-J 

Ui 

it 

X 

X 

X 

X 


X 



Ui 


X X X Z} X 


4-4 

a -o 

\/\ 

4-4 

4 — » 


G 

CO 

«’/) 

G X 

X 

! x 

J 

M 

X 


... 

•-■ 



'w 



p 


*w* ••-•' ^ w 

•*-* 

+ 

1 — ’W 


+ 

t— 

©• 

H- 

4—4 


H* ’»■ 



** 


21 











iri 


X 



X 

X 

11 

H* 





2 

2 

u 

i— . 

X 

X 

X 

X 


X 



IJ 


LlU.Ll 1 CU.LL 

11 

C X X 

II 

a 

lyj 

X 

a x 

I u. 

X 

1 1 



X 


«— 4 

t — 4 

1—4 


1—4 



V -4 


i —4 *— » tV'i *X 

1—4 

X 

in 4 -t 

t— 

X 

t’X 


1 — 

© 

t —4 

h— t * K ^ 


£ 

1/1 

in 







•si* 


I£| 















•I 

*"• 

□ 










a 















U; 

Xr = 

y 










J 


rf 00 



o 



OJ 





OJ 


i.M 

1 - 

ij 

4—4 

I 1 








1 f 

IJ 

1 1 

CU OJ 



0 J 










IJ 1 1 

! 


IJ 

-b-4 CU 

rr, 

rr 

uo 

X 

r- 

co 

iT. 

IJ 

IJ 

OJ 

co rr uo x r- 

CO 

tT. 

■r 4—4 

OJ 

co 

rr 

uo 

X' 

r-- 

CO 

iT, ©• 

4- 

Cu C ri 


2 - 

'IT- i X‘ 


>”• 

*“ 


o 


o. 


4 — 

4-4 

—H 4-4 4-4 4-4 4 — 

—4 

4-4 

OJ cu 

cu 

OJ 

cu 

'OJ 

OJ 

cu 

OJ 

OJ CO 

co 

0 *' 0 " 


5 * 

*r nr 

Tt 

rr 

rr 

rr 

rr 

rr 


rr 

rr 

rr 

TTTTT 

rr 


rr rr 

rf 

rr 

rr 

rr 


rr 

rr 


rr 

’ T t; 

rr 


; 2 ' •— 1 



2 ' 

•2 

2 ' 


o 


• 2 : 

2 ' 


2 - 


2 ' © 

2 * 





2 > 

© 

©' ©• 


2 ‘ © 

■2 


96 



Eb <K> . bT . THTM I H ( I > "> THTM I H < I > =SEb <K> 


t — i iT» CO 


X X 


H- 

3 

E 


UJ 

E 

X 

Ml 

UI 

E 

Ml 

Ml 

UJ 

E 


Mi 

UI 

J Mi 
X 2 
1. □ 

ft ft* 


J K 
X Q 

p 

X Ml 

Z « 

lL > 

□ y 

Ui c 
M‘ ft 
UI 2 

lj it 

□ ft 
E UI 
Cl E 








E 


« 

Ui 

x 

UJ 

Ml 









Mi 

ft 

t t 

X 

H® 







ij 


E 

x 

a 

ft 

X 







M 


Ui 

j 

j 


ft 







ft 


ft 

a 


X 








x 


Ui 

Cl 

c 

in 

$ 







i- 


z 

E 

ft 









Ml 


c 

Ui 

Ui 

□ 

E 









E 

ft 

X 

E 

x 







a 

*Z* 

X 

z 

ft 

X 

Ml 








• 

E 

HO 


ft 

Mi 








<hH 



IX 


Ui 

» 






X 

1 


z 

ft 

Q 

IJ 

UJ 






UJ 

X 

D 

H* 

111 

in 

Ui 

E 






X 

CO 

H» 

€ 

E 


z 







4fc 

cu 

ft 

ft 

u 

2 


Ml 







• 

E 

ft 

Ml 

2 

r* 

Ml 






E 

♦ 

□ 

D 

HO 

HO 

a 

UI 






UJ 

♦ 

ft 


M 

z 

z 

E 






ft 


Ml 

a 




E 





L 

z 

.H., 

HO 

ft 

ft 


Mi 






X 1 

ijn 

ft 


c 

E 

H* 

u 






X. 

> 


iu 


E 


HO 







X 

J 

z 

in 

Ui 

ft 

ft 

or* 





I 

X 

c 

H* 

Z 

ft 

HO 

c 

nj 


CO 


Z* 

IJ 

CO 

M 

ft 

HO 

z 


ft 




• 

x 

\ 

ft 

z 

E 

HO 

I 

Ml 

a 


a 


hH 


in 

x 

Q- 



(£( 


ft- 


ft- 

HH 

II 



E 

E 

X 


ii 1 

Ui 


HH 




UJ 

x x 


ft 

i i 

Ui 

□ 

in 

□ 


□ 

LL 

HH 

in 

X '2 

Ui 

Z 

i 

z 

X 

<1 

IJj 

<*H 

CD 

X 

A 

»t 

ft- x 

J 

h't 

Ui 

HO 

ft 

E 



Cl 


LL > 

E 

2Z 

Cl 



ft 


UJ 


Cl 


LL 

X hh 

'I 

co x 

T* 

1 z 

z 

2 

z 


ij 

uo *x r- r 

ft ft ft 








HH 


+ 

*w 

a ft 



A 


% 

A 



M 



CO 

ft- 


a 

D 

r- 

u» 

2 


® 


£ 



ft- 


X 

ft- 

ft- X 

ft- 

• 

\ 


£ 

X 

iVi 


Ml 


A 


X 


M 


E 

2 


• 





X X 

CO 

X X X CO 

X 

“j 

X 




A 




X 

H 

ft- 


H 

M 

m 

2 

CO 

X' 


lL 



ft- 

2Z 

ft C ft Z ft □ 

9 

co 


\L 


X 


E 



11 



X 

ft- 

Z 



11 


□ 



o i> ;« 


X 

1 

CO X 

i 

q; 

ft- 


□ 

ft 



□ 


H 

ft- 

HH 


2 

Z 

Mi 

z 

ft- 






■ a 

f—o 

ft- X 

• ^ ft 

LL X 

X 



•X 

ft- 




11 

X 

4^ 


C 

ui 


a 

X 

>- UJ 


21 



® 1- 

11 

CO 

ft 

ft- il CO 

ft 

• 

ft- 


z 

• 

\JL 


Z 


X 

ft- 

o uo 


E 

z 

ft 

HS 

Cl 



□ 


Z r 

11 X 

HH 

11 


a -> u 

— »m 


cO 


□ 

uo 

•Z‘ 


□ 


PH 

HH 

a ft 


E 

UI 

X 

ft 

X 

T" T* 


HO 


• 



ft 

HH 

X X 

HH 

HH 

11 


H4 

fl 

iV‘» 


HP 




A il 


UJ 

E 

X 

X 


pH 


ft 


Z‘ 

r/» • 

o 

(A 

A 

- CU CO 

Lh 

A 

t£i 


ft 

X 

11 


ft 


z 

Z : 

H— * iJL. 


ft 

2 

ft 

2 

-J 

X ft* 


X 


II 


( V( 

ft- 

H 

CO ft- 

11 


“"2- 


*1 

O 



c 


ft 

ft 

X ft 



hi 


G 

-J 

i— X 


2 



X 


X 


X 

— j. 


X 


2 

X 

X 


2 




X X 


t-H 

•r 

Ui 

E 

x x a 


ft 


HH 

ft- x a 


HH 

X C ft- 

HH 

X 

X 


Z* 

X 

X 


z 


X 

a 

Il X 



E 

ij 

r> 

Uta, 

i2 

ft- <2 


IJ 


A 

(/} M 

A 

co 

A 

»— * L—0 CO 

A 

HH 

co 


1 1 

X 

X 


n 


1—4 

A 

X ft 


2 

Z 

H. 

2 




2 












2 




2 






2 


ft* 

ft* 




•I 












»r 




C 






€ 

E 

Q 

Z 




Ij 





t — , 


iT’i cu 





ij 




Ij 






lj 

a 

Z 

H. 

00 

<T, 

u 

IJ 

IJ 




0*» 


CU CO 


CO 


IJ 

ij ij 



IJ 

lj 

IJ 




Ij 

lj 

u 

ij 

i J 

•X 

— * CU 

|V. 

ft 

UO 

'X 

r - ‘- CO 

CT\ 

(~i 


CU CO TT 

uo 

•L 

r*. 

co 

)T> O 

*—4 

cu 

co 

ft 

uo 

*£• 

r-. 

CO fft 

Z- 

HH 

Cu 

C r ' 

ft 

ft 

ft ft 

ft 1 

ft 

ft- ft- 

ft- ft- 

ft- 

uo 

UO 

uo uo uo uo 

uo 

uo 

uo 

UO 'X* 

l X 

‘■X 

•X* 

•X' 

•X 

•X 

*x« 

■X X 

r- 

r-- 

r- 

r*-- 

r-. 

X 

ft 1* ft 

ft 

ft- 

ft- ft- ft- 

ft" 

ft 

ft 

’i- rr ft- 

ft- 


ft* ft- ft- ft- 

ft* 

ft* 

ft 

ft 

ft 

ft 

ft 

ft ft 

ft 

ft 

ft 

ft 

ft 


*—■ w 

*— • 

w 

Z 


' w 

® 

o 

Z‘ 

O 1 


•z* 



’ w 

*— 1 1 




o 

Z‘ 







■— 


97 



0475c EXTREME SINCE THE RBITmL I NTERPDLAT I DM IS NOT EXECUTED TOO MANY TIMES 

0476c 

0477 CALL IHTERP 

0478 4 0 RAPP < I THT , I R> = RVALUE 

0479c 

0480c CALL- SUE ROUT I ME TO CALCULATE SPECIFIC DISTORTION PARAMETERS. 

0481c 

0488 CALL D I SPAR 

0483c 

0484 RETURN 

0485 END 

0486c 

0487c 


UJ 

CD 



0488c 

0489 SUBROUTINE EKTRME <S6K > KSS » KMAX» T » FLUR> 

0490c 

0491c SUBHDUT X HE WHICH MANAGES COMPUTATION OF EXTREME VALUES 

0492c S5KF IS DISTORTION FACTOR SIGMA 

0493c KSS IS MEAN VALUE 

0494c KM AX IS EXTREME VALUE 

0495c T IS TIME ON PD I NT » SECONDS 

0496c FLUR IS VORTEX FLUX RATE. 

0497c 

0498 REAL. KSS»KMAX 

0499 S6KF=S6K 

0500 110 fi K S = K S S / S 6 K F 

0501 FLUX=T+FL.UR 

0508 DKS=0. 0 

0503 1 HBX=0. 0 

0504 IF (SGKF.EQ. 0. O') 60 TD 120 

0505 IF <BKS.6T.64. 0> BKS=64. 0 

0506 CALL SEARCH <AKS> DKS» FLUX> 

0507 K M B X = K S S + D K S ♦ S 6 K F 
05 08c 

0509 120 RETURN 

0510 END 
0511c 

0512c 




0516C SUBROUTINE TO TRANSFORM INPUT MEASUREMENTS INTO MAPPING PARAMETERS 

0517c SUBROUTINE ALSO CONTROLS SOME OUTPUT 
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1638 92 CALL F'F'HT <£> 

1639 CRLL F'FIX < . TRUE 

1640 CALL DISTRT 
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HIS SUBROUTINE COMPUTES THE RFiDIhL. AND CIRCUMFERENTIAL VELOCITY 
RhI'IENTS mT EmCH grid point using cubic spline interpolations. 
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THIS SUBROUTINE PERFORMS CUBIC SPLINE INTERF'OLhT IONS (SLDFE“BmSED) 
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APPENDIX B. 


OPTIONAL "SEGMENTED VORTEX” ADDITION 

An addition to the source code given in Appendix B 
includes the Segmented Vortex approach described in Ref . 8 . 

In this approach, the vortex model derived in the standard 
Melick approach is divided into eight segments. This process 
allows for simulation of a nonlinear vortex, or a vortex 
ring. Theoretically, this should allow for more accurate 
modelling of the inlet flow distortion, for better results,, 
The segmented vortex approach showed some improvement of 
the predicted peak instantaneous distortion contour map in 
certain cases [Reference Q) . The User is free to experiment 
with this addition to the source code. The following page 
lists the addition, showing where it is to be inserted in 
the original code. 


149 



IF (RAD.GT.o70.u) RAD=670.0 


: DP T 1 IS THE PRESSURE FLUCTUATION AT PPQBE A R ( J ) AND ANGLUC(J,1 

DPT1=2.0*V9aP*BaSE/£XP(0.5*KAd) 

75 P(J,I)=HS(J/I)+L)PT1*QUPT*TMAVGG 
80 CONTINUE 

00 TO 90 . . , 

original 

RESET STATIC PRESSURE kTNGS TO STATIC PRESS"RF. code 

'78 00 7R 1=1, NR 4 

79 P(J, I)=PS C J/ n 
90 CUNiI.MUE 

*** NEW "SEGMENTED" CORTEX METHOD 8Y STFPH'EN DFnNOW *~*7 I 
FIND PROBE LOCATION iRaDLOC) NEAREST TO MF AN LINE 

OU 120 I=1,NP 
0 0 P ( I J = 10. 

DO 1 10 J=i,oR 

DP ( J , I ) = AbSCTMAVGG - PSCJ/in 
IF ( D p ( J , I ) .LE. DDPCI)) GO TO 100 
GU TO 110 

ioo oup(i) = jPCJ,n 

RADLOh(I) = RADLOC(J) 

110 CUNTINUE 
120 CONTINUE - 


SOLVE FOP JEN PEAK INSTANTANEOUS PRFsSUPE APRAY 

DU 190 1=1 , iP 
00 150 J = 1 / N R 

TF (J .lT. J1 .uR. J .GT. J2) GU TO 1 25 
DRR = A d S (RaDlOMC I) - RADLOC(J)) 

PR A = DRR/R T/AR T 
PAD = RpA**2-l. 

DRTi = ARS C2*V«AR*PKA*QOPT*T M AVGG/EXP( ,5*R AU) ) 
IF (PS(U,n .LE. T M A V G G ) GO TO 121 
P(J,I) = PS(J,T) + DP T 1 


addition 


GQ TO 130 

121 P(J, I) = PSCJ, I) - UP T 1 
GO TO 130 

125 PI J# I ) = PS( J r I ) 

130 CONTINUE 
190 CONTINUE 

IF ( I V - 1 ) 91,9 1,92 

91 CaLL pfix C.TkUE.J 

call [jIstrt 

GU TO Boo 

92 CALL PRNT (2) 

CALL PFIX C . TRuE. ) 

CALL OISTRT 

CALL PFIX (. FALSE. J 
CALL rlAINLP 
CALL LNpOUTtl) 

WRITE (o,1995) 

WRITE (6,2000) vraR,art,g,b,ayrt,vl,rrt,tv 

WRITE (o, 1999) VBM/VBO/AO 


oniginal 

code 
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APPENDIX C. 


QUICK REFERENCE GUIDE 


MAXIDYN Program Description 
Input Data Description 
Output Data Description 



1 . MAXIDYN Program Description 

The MAXIDYN dynamic distortion program computes the 
most probable peak instantaneous distortion level given the 
steady-state distortion conditions, and generates a peak 
distortion map based on the predictions. The Melick corvee- 
ting vortex model and statistical approach is used in this 
predictive analysis, with some modifications and improvements 
to enhance program Flexibility. The complete FORTRAN program 
requires sufficient computer memory capacity for approxi- 
mately 100,000 words, plus typically 5000 words per data 
set. Run time varies from system, but is generally limited 
only by the online printer output capacity on most main- 
frame systems. 

Input data requirements include the rake and probe 
conf iguration used in the test, the steady-state static 
and total (stagnation] pressure measurements in the rake 
plane, some basic inlet flow parameters, and optionally the 
root mean square pressure fluctuation level measurements. 

The input data are described briefly in the next section, 
and formatting requirements are shown in the figure. Some 
of the input data are optional, that is they may be deleted 
from the input file. The program automatical ly assigns pre- 
defined default values, or as in the case of the rms levels, 
the data are computed based on other input data. 

The printed output of the MAXIDYN program includes 
several pages of steady-state and dynamic distortion data. 

The input data are organized in groups printed in the first 
few pages, along with steady-state distortion computations, 
Melick vortex model parameters, and the steady-state map. 

The rms fluctuation levels and/or turbulence modelling data 
are printed in the next few pages, along with the statistical 
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predictions of the most probable peak instantaneous distor- 
tion levels. Finally, a dynamic distortion map is generated 
based on the peak instantaneous prediction. A brief descrip- 
tion of the output data is given in the third part of this 
Appendix . 

2 . Input Data Description 

The table below briefly defines the variables in the 
input data. Certain data may be deleted from the input file 
because they are considered optional and generally are not 
part of the computational procedure . Some of these data are 
assigned default values as needed within the program. All 
optional data are indicated with a in the second column 

of the table below. In the Format column, Tt F TT indicates a 
real number, TT I TT indicates an integer, and TT A Tt indicates an 
alphanumeric array, according to standard FORTRAN rules. The 
arrangement of the input data is illustrated in Figure A1. 

Data * Format Definition 

NR 15 Number of pressure tap radial locations 

NP 15 Number of rakes 

RAOLOC F10.5 Radial location of pressure taps 

ANGLDC F10.5 Angular location of rakes in degrees 

KD 15 Distortion factor selection index 

KD = 1 : KTHETA 
2: KD2 

3: IDC 

4: IDR 

5: KRA 

6: KA2 
7: DSPR 


8: ID 
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Data 


Format 

Definition 

T 

-j* 

F10.5 

RMS fluctuation measurement time-on-point 




default = 1 sec. 

FO 

* 

F10.5 

Engine filter frequency, default = 500 Hz 

FCO 

* 

F 1 0 . 5 

RMS cut-off filter frequency, 1000 Hz def 

TITLE 


A80 

Title block 

BRP 

T" 

F10.5 

Base radial profile array, default = 1.0 

PS 


F10.5 

Steady-state pressure array 

ALPH 


F5 . 3 

Inlet angle of attack in degrees 

PSI 


F5 . 3 

Inlet yaw angle [crasswind] in degrees 

MO 


F5 . 3 

Freestream Mach number 

UP 


F5.3 

Inlet flow velocity in fps 

BF 

* 

F5.3 

B-f actor for weighting in KA2 computation 




default = 1.0 

CKP 

* 

F5. 3 

Circumferential ID weighting factor. 




default = 16.4 

RKF 

* 

F5.3 

Radial ID weighting factor, default = 11. 

MFR 


F5. 3 

Inlet mass flow ratio, default = 1.0 

NTUR 


15 

Dynamic data selection index 

NTUR = 0: Dynamic data ( rms levels) input 




1 : Dynamic data synthesized 

SPTRC 


F5.3 

Supersonic inlet pressure recovery. 




default = 1.0 

The following data 

is required if NTUR = 0: 

NPG 


15 

Data run number 

NPR 


15 

Dynamic probe location number 

RS 

* 

F5. 3 

Filtered rms level -J* unfiltered rms level 




default = 0.5 

SIG 


FB.4 

Unfiltered rms fluctuation level 
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rORTRAN STATEMENT 
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Input: Data File Formatting Arrangement 
























3. Output Data Description 


The printed output of the MAXIDYN program consists of 
several pages of data and computations. The first two pages 
are related to the steady-state distortion and some Melick 
vortex model parameters, the middle page or pages are 
related to the dynamic data and the statistical determination 
of the most probable peak instantaneous distortion, and the 
last two pages are related to the generation of the peak 
instantaneous distortion map. The contents of each of the 
pages of output are briefly defined below. 

Page 1 


At the top of the first page, immediately below the 
title blacks supplied by the user, is the steady-state 
PRESSURE ARRAY. This array is identical to the pressure 
array of the input data, except the static pressures have 
been deleted. Each column of the array represents a rake , 
while the rows represent probe locations of probes along 
the rake. Immediately below the pressure array is the BASE 
RADIAL PROFILE, also from the input file. The next table, 
the OVERALL FLOW DESCRIPTORS, provides some simple distor- 
tion parameters: 

PTMIN is the minimum measured local pressure from the 
steady-state PRESSURE ARRAY. 

PTMAX is the maximum pressure from the PRESSURE ARRAY, 

PTAVG is the average pressure from the PRESSURE ARRAY. 

PSAVG is the average static pressure . 

QAVG is the average dynamic pressure. 

The next table is a listing of the eight FLOW DISTORTION 
FACTORS and their values, along with weighting factors, 
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At; the bottom of the first page, some of the Melick VORTEX 
PROPERTIES are given. These properties are: 

THMN - The rake showing minimum average pressure 

RKMN - The average pressure along rake THMN 

THMX - The rake with maximum average pressure 

RKMX - The average pressure along rake THMX 

□TH - The angular difference between THMX and THMN 

THETA - The angular location of the center of arc DTH 

A 1 - The radius of the steady-state Melick vortex 

G1 - The orientation angle of the steady-state vortex 

Page 2 


The next page of printed output is a distortion contour 
map for the steady-state case. Relatively high and low pres- 
sure regions are indicated by symbols, which represent the 
percent difference from the average pressure, as indicated 
by the KEY TO MAPPING SYMBOLS immediately above the map. 

The average pressure is printed to the right of the map. 

Page 3 

Page three of the output includes a listing of some 
of the input data, including flow parameters, and the 
dynamic data, assuming dynamic data was included in the 
input file. If dynamic data was not included in the input 
file, the third page would include some data from internal 
turbulence calculations. The data on this page includes: 

T - The dynamic data time-on-point, in seconds 
FO - The engine filter frequency, in Hz. 

RT - The outer rake diameter at the static tap 
RI - The centerbody hub radius 

SPTRC - The supersonic inlet shock pressure recovery 
ALPH - The inlet angle of attack 
PSI - The inlet yaw/crosswind angle 
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MO - The freestream Mach number 

ETA - The average pressure at; the measurement: plane 
MFR - The inlet mass Flow ratio 

U2 - The inlet Flow velocity at the measurement plane 
QPT2 - The ratio oF the dynamic to total pressure 
FC - The rms Filter cutoFF Frequency 

IF dynamic data is included in the input data File, these 
data are printed in a table. The terms in this table are : 
PROBE - The dynamic probe location index 
RS - The ratio aF Filtered-to-unFiltered rms level 
SG/PT2 - The unFiltered rms Fluctuation level 
A/RT - The vortex radius resulting From the rms level 

The rms Fluctuation level and vortex size are given below 
the dynamic data table. 

IF dynamic data are excluded From the input File, these 
data must be synthesized by the turbulence modelling scheme. 
In this case, the dynamic data table is replaced by a table 
oF DIMENSIONLESS VELOCITIES occurring at each oF the probe 
locations, and a table oF iterations oF turbulent kinetic 
energies. These are provided For the convenience oF the user 
and are not directly involved in the distortion analysis. 

Page 4 


The Fourth page oF output includes a listing oF the 
DISTORTION FACTOR EXTREME VALUE computations. For each of 
the eight distortion Factors, values oF the most probable 
peak instantaneous distortion are presented. This table 
includes the STEADY STATE, MEAN INSTANTANEOUS, and peak 
instantaneous distortion at various statistical conFidence 
levels. In addition, the Filtered and unFiltered rms dis- 
tortion levels are indicated; INF reFerring to the unFiltered 
case and FO representing the Filtered case. 
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In *the case where the turbulence modelling scheme is used to 
generate the dynamic data, the fourth page includes further 
results of the TURBULENCE CALCULATIONS. This table is similar 
to the dynamic data table as described in the third page of 
output, with the exception that the term RS is deleted, and 
velocity gradiants ( UU ) and turbulent kinetic energy terms 
(K and E] are added. The unfiltered rms levels are presented 
in the last column. The DISTORTION FACTOR EXTREME VALUE 
table is moved to the fifth page in this case. 

Page 5 

The fifth page is arranged exactly like the first page, but 
with notable differences. All of the terms in the PRESSURE 
ARRAY, OVERALL FLOW DESCRIPTORS, and FLOW DISTORTION FACTORS 
tables refer to the peak instantaneous case rather than the 
steady-state case. In addition, the VORTEX properties table 
contains additional terms: 

VBAR - The vortex "strength”, or maximum swirling 
velocity 

A/RT - The radius of the Melick vortex 
GAMMA - One of the vortex orientation angles 
BETA - The second vortex orientation angle 
AY/RT - (The same as A/RT] 

VL - The vortex length limit (generally "infin ity"} 

R/RT - The radial location of the vortex core 
VBMX - The maximum instantaneous vortex "strength” 

VBO - The approximated mean vortex strength found in 

an iteration of strengths and distortion factors. 
AO/RT - The vortex size indicated from the rms data. 

Page 6 


The sixth page is similar to the second page except the 
distortion map is for the peak instantaneous case. 
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VI. CONCLUSIONS S RECOMMENDATIONS 


The subject computer program can be used to aid the 
prediction of maximum instantaneous distortion levels, and 
the peak instantaneous contour map, given the steady™ state 
distortion data and, optionally, the dynamic rms pressure 
fluctuation data. There are some improvements which can 
be added to the program, at User’s discretion. 

One improvement currently being researched at the 
University of Kansas is the replacement of the single 
steady-state vortex model with a series of vortices whose 
axes lie approximately along the ’’mean line TT of pressure 
recoveries at the compressor face. This effort is intended 
to improve the predicted peak distortion contour map to more 
closely resemble the experimental map produced by DYNADEC . 
Other methods of improving the accuracy of both the peak 
distortion level, and the corresponding contour map, with 
respect to experimental results, would be highly desireab le . 

The accuracy of the present analysis is discussed in 
References 5 and 6, along with the basic derivations in the 
theoretical analysis. In general, the Melick technique is 
reasonably accurate for preliminary design and analysis „ The 
major benefit of the Melick method is it’s efficiency, and 
the general tendency to over-estimate the experimental or 
true peak distortion level, rather than under-estimate it. 
One of the primary difficulties with the Melick analysis is 
in predicting the distortion levels for inlets with separa- 
ted flows. It would be desireable to try to improve the 
accuracy of the peak prediction for this extreme case, which 
can occur especially often in highly maneuverable aircraft, 
which operate at high angles of attack and yaw angles, and 
also tend to have complicated inlet duct shapes. 
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